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ABSTRACT 


The Russian designed and manufactured TOPAZ-II Thermionic Nuclear Space Reactor has 
been supplied to the Ballistic Missile Defense Organization for study as part of the TOPAZ 
International Program. A Preliminary Nuclear Safety Assessment investigated the readiness 
to use the TOPAZ-II in support of a Nuclear Electric Propulsion Space Test Mission 
(NEPSTP). Among the anticipated system modifications required for launching the 
TOPAZ-II system within safety goals is for a U.S. designed Automatic Control System. The 
requirements and desired features of such a control system are developed based upon U.S. 
safety standards. System theory and design are presented in order to establish the basis for 
development of a hybrid control model from available simulations. The model is verified 
and then used in exploration of various control schemes and casualty analysis, providing 


groundwork for future Automatic Control System design. 
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I. INTRODUCTION 


A key element among the long term technological requirements for advanced space 
missions ıs the required development of a suitable high power electrical generation system. 
A joint National Aeronautics and Space Administration (NASA)/Department of Defense 
(DoD) review committee determined that Space Nuclear Power (SNP) is the best option 
for meeting demands for certain high power and long hfe applicanons. The Strategic 
Defense Initiatve Organization (SDIO), now the Ballistic Missile Defense Organization 
(BMDO), led the effort to acquire Soviet technology in order to advance SNP science 
within the United States. Within this process, the Russian TOPAZ-II space nuclear power 
system has been part of a joint testing program since December 1991. (Malloy et al., 1994) 


The TOPAZ-II system provides power using thermionic energy conversion. This 
process is a form of direct conversion in which heat is converted into electricity without 
intermediate energy transformation stages. Unlike conventional heat-power systems, 
thermionic conversion requires no moving mechanical parts. This eliminates the need for 
frequent maintenance, making it particularly suited for use as an electrical-power source for 
deep-space probes and various other space vehicles. Thermionic systems have low system 
weight and the ability to withstand acceleration higher than 30 times the force of gravity. 
The heat for thermionic devices may come from any convenient source. The TOPAZ-II 
design uses nuclear fission heat generation, which allows a high power level over long 


system life, independent of location in the solar system. 


The SDIO initially intended use of the TOPAZ-II reactor for a Nuclear Electric 
Propulsion Space Test Program (NEPSTP) Mission but launch of the TOPAZ-II requires 
determination that the mission could be performed safely and with reasonable costs. A 
Prelimmary Nuclear Safety Assessment performed by a team organized by SDIO 
determined that modifications would be required to meet U.S. safety requirements. Among 
these changes 1s the requirement to design and build an Automatic Control System (ACS) to 
U.S. standards. The Russian ACS 1s a relatively massive analog system that cannot be flight 


qualified and does not incorporate certain desired features. (Flight Safety Team, 1992) 


Additionally, the TOPAZ-II system could be made more robust through minor 
alteration of system design. Certain mission profiles for space nuclear power applications 
include operation for long periods at reduced power levels. Recent research (Benke, 
Venable, 1995) demonstrated work on the road to feasibility of the TOPAZ-II system for 


use in such missions. 


The goal of this work 1s the initial design of a start-up control system that meets all 
safety requirements and is robust enough to enable the expanded mission possibilities of the 


TOPAZ-II reactor. 


II. BACKGROUND 


The successful exploitanon and exploration of space will require ever increasing 
amounts of electrical power. Future missions that will demand large power sources include 
high-bandwidth communications and remote sensing satellites, deep space exploration, 
extraterrestrial manufacturing, manned missions to Mars, and colonization of the Moon. 
Many of these tasks will take place in harsh. environments where solar power is limited, 
intermittent, or unavailable. Nuclear power and thermionic conversion can serve as a 
compact, durable energy source that is capable of long-term operation in radiation belts, 


deep-space vacuum, and hostile atmospheres. 


Direct energy conversion is a process by which heat is changed into electrical energy 
without the use of intermediate stages or moving mechanical parts as required in such 
conventional systems as the Rankine cycle. The two main methods of direct conversion 
incorporated with a nuclear heat source are thermoelectric and thermionic. In a 
thermoelectric power generator, a temperature differential between the upper and lower 


semiconductor material being generates power. 


The thermionic method uses two electrodes. The emitter is heated to a sufficiently 
high temperature to emit electrons. The collector, which receives the emitted electrons, 
functions at a much lower temperature. Thermionic converters can operate at significantly 
higher radiator temperatures than the thermoelectric process; resulting in smaller radiator 
size and mass. This higher temperature range potentially leads to much higher overall 
efficiency. Thermionic conversion systems also restrict the high-temperature region to the 
fuel and the emitter. This permits the use of lighter and cheaper low temperature structural 
materials. (Paramonov and El-Genk, 1994) 


Both the United States and the former Soviet Union have conducted advanced 
research on direct energy conversion systems. In 1965, the U.S. launched its only 
experimental thermoelectric nuclear reactor space system ever flight-tested, the SNAP-10A. 


The Soviet Union had an extensive space nuclear power program that included 35 nuclear 


reactors; incorporating both thermoelectric and thermionic methods. TOPAZ-II ıs the 
most technologically advanced single-cell thermionic system ever built. By March 1994, all 
TOPAZ-II space power reactors in the Russian inventory have been brought into the 
United States for use in the TOPAZ International Program (TIP). Researchers from Russia, 
France, the United Kingdom, and the United States conduct nonnuclear testing at the New 


Mexico Engineering Research Institute (NMERI) in Albuquerque, New Mexico. 


A. SPACE NUCLEAR POWER APPLICABILITY . 


The dominant factors that determine the choice of electrical power source are the 
electric power level required and planned use duration. Nuclear power 1s particularly suitable 
for applications that demand high power levels over a long period of time. Typical examples 
of such missions are for operation of large space platforms, lunar bases, high power 
communication and remote sensing satellites, specialized deep-space probes, and a Mars 
landing. Thermionic systems with nuclear reactor heat sources are optimal for electrical 


requirements in the range of 10 kW to 1000 kW (Larson and Wertz, 1992). 


Nuclear power has several advantages over conventional power. These systems 
typically have a compact size with a low to moderate mass. Nuclear power may be provided 
over a long mission lifetime in hostile environments, independent of the sun. A nuclear 
reactor system exhibits very high hardness with respect to natural nuclear radiation, nuclear 
or laser threats, and meteorite or pellet impingement. These systems exhibit high stability 
and maneuverability. However, there are many concerns with respect to the use of these 
systems. Chief among these considerations 1s the potential environmental impact, which 
imposes the requirement for a strict safety analysis reporting program (Larson and Wertz, 
1992). Additionally, space reactor and direct energy conversion technology is less developed 
than either chemical or solar power. Past nuclear missions have included several important 
space programs including Apollo and all major interplanetary probes (Benke, 1995). In sum, 
space nuclear power 1s a vital niche technology whose importance will grow with our future 


exploitation of the solar system. 


B. SPACE NUCLEAR POWER PROGRAMS 


Space nuclear power has been under development since the birth of space 
exploration. The majority these programs have been conducted in the United States and the 


former Soviet Union. 


1. United States Space Nuclear Power Programs 


A favorable political environment and high public support for nuclear power 
technologies fostered its development from the mid-1950’s to the early 1970’s. The first 
major space use of nuclear power was in the Transit 4A navigational satellite, launched in 
June of 1961, which employed the Systems for Nuclear Auxiliary Power (SNAP) version 3B 
source. This power system relied upon a Radioisotope Thermoelectric Generator (RTG), a 
method that uses radioactive decay processes as a heat source. RTG’s are highly reliable as 
proved in the performance record of 24 U.S. missions, including Pioneer, Viking, and 


Voyager (Benke, 1995) 


The 1955 Project Rover program at Los Alamos National began research of 
space nuclear technology. This program was the start of 18 years of nuclear rocket research 
and development that advanced the search for a suitable high temperature reactor. The 
Nuclear Engine for Rocket Vehicle Applicaton (NERVA) program demonstrated the high 
thrust and specific impulse capabilities of a high power density nuclear rocket (Angelo and 
Buden, 1985). 


U.S. development of nuclear reactors for space electric power applications 
began in 1957 with the production of the SNAP-2 reactor. This system combined a thermal 
uranium zirconium-hydride reactor as a heat source with a liquid mercury Rankine cycle for 
the thermal to electric power conversion. Following the SNAP-2 system was parallel 
development of SNAP-8 and SNAP-10 programs. The SNAP-8 system was essentally a 
scaled-up version of the SNAP-2 system, providing a higher power level. The SNAP-10A 


system used the thermoelectric conversion method in place of the Rankine cycle. 


In April 1965, as part of the SNAPSHOT mission, SNAP-10A was the first 


reactor ever launched and orbited, and the last such system flown by the U.S. The 


SNAP-10A system generated approximately 600 watts of electrical power in an orbital 
alutude of about 1300 km. This reactor operated at power for 43 days before the on-board 
Reactor Control Unit (RCU) safety system initiated an automatic shutdown. An apparently 
false detection of spacecraft voltage regulation failure caused this event. The downlinked 
spacecraft data indicated that system operation was normal prior to this event. (Malloy, 
1994). This incident highlights the importance of RCU design to overall system reliability 
and safety. 


After major advances in solar power, space nuclear reactor technology 
languished in the 1970’s. Thermionic research had developed to the point that it had 
supplanted thermoelectric development by 1969. American researchers tested several 
thermionic devices; however, no complete reactor system using thermionics has ever been 
built ın the United States. The U.S. space program continued to use nuclear power sources, 


mainly plutonium fueled RTG's. 


During the 1980’s, the inauguration of the Strategic Defense Initiative (SDI) 
revitalized space nuclear research. The SDI program spurred the development of the Space 
Power Advance Reactor (SPAR), later renamed SP-100, to meet projected 100 kW electric 
needs. The Department of Energy (DoE), NASA, the Defense Advance Research Projects 
Agency (DARPA) established a combined program. These organizations conducted trade- 
off studies comparing thermoelectric, thermionic, Rankine, Brayton, and Stirling energy 
conversion methods. The lack of reliability and performance data for thermionic systems 
within the U.S. resulted in the choice of a thermoelectric design for the SP-100 system. In 
order to exploit the superior power to mass ratio of in-core thermionic power conversion, 
DoD and SDIO began the Thermionic Fuel Element Venfication Program (IFEVP) in 
1985. Development of the SP-100 system has essentially been discontinued with the change 
in the political situation following the end of the Cold War. (Malloy et al., 1994). Figure 1, 


below, shows a drawing of the SP-100 mission concept. 
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Figure 1. SP-100 U.S. Space Reactor. 


2 Former Soviet Union Nuclear Power Programs 


Soviet space nuclear power design began in the 1960's, with several systems 
that were largely derivative of early U.S. efforts (Malloy et al, 1994). The Soviet space 
nuclear power development benefited from a more constant level of effort, resulting in a 


total of 38 reactors launched over the history of its space program (Benke 1995). 


The first Soviet space nuclear power system, the BOUK reactor, was 
launched in 1967. This reactor configuraton utilized thermoelectric conversion in a 
configuration similar to the U.S. SNAP-10A system. This 1-2 kWe reactor design was 
extensively used in the Radar Ocean Reconnaissance Satellite (RORSAT) program, which 
launched an estimated 33 satellites. Nuclear power was the preferable source since a solar 


powered RORSAT would have limited to an orbital lifetime of a few days due to the effects 





of atmospheric drag on the large solar panels. The orbital life of the RORSAT system was 
approximately 100 days. (Malloy et al., 1994). 


The inadvertent reentry of COSMOS-954, a RORSAT system, caused the 
most well known space nuclear power incident. In 1977, this satellite dispersed reactor 
material over a 600 km path in northern Canada (Angelo and Buden, 1985). The end-of-life 
orbital boost system did not activate due to a telemetry fault (Malloy, 1994). Following this 
mishap, the United Nations adopted international nuclear space reactor design and 
operating principles. In order to achieve a “nuclear safe” orbit, improved RORSAT design 
automatically separated the reactor core from the satellite and placed it in a higher altitude 
of approximately 950 km at the end of mission life. (Further discussion of nuclear safety 


considerations 1s found in Chapter V, Safety and Functional Requirements.) 


The Soviet drive for higher altitude surveillance lead to the development of 
the COSMOS systems. The two reactors, COSMOS-1818 and COSMOS-1867, were 
launched in flight tests conducted in February and July of 1987. These satellites were test 
platforms for an advanced reactor design utilizing mult-cell thermionic conversion 
techniques. Both tests were extremely successful, lasting for 142 days and 342 days 


respectively (Paramonov and El-Genk, 1994). 


A mult-cell thermionic fuel element (TFE) uses a stacked series of short 
thermionic cells, rather than the long unitary fuel element of the single-cell design. The 
TOPAZ Space Nuclear Power System (SNPS) (translated from the Russian designation: 
TOJIA3 Kocmnyeckaa flaepiaa Sueprernyeckas Ycranoska — KA3Y) developed multi- 
cell technology. The ENISY program produced a single-cell design. American researchers 
incorrectly designated the ENISY program as TOPAZ-II (Voss, 1994). TOPAZ itself 
subsequently became known as TOPAZ-I in the U.S. No single-cell TFE system was ever 
flight tested, although 26 ENISY reactor systems were built. Testing of 19 of these systems 
took place between 1970 and 1989 (Schmidt et al., 1994). 


C. TOPAZ INTERNATIONAL PROGRAM 


The TOPAZ International Program (TIP) was an outgrowth of the reduced 
tensions at the end of Cold War. Both the U.S. and Russia faced large funding reductions 
and program cancellation in all space related nuclear research. With the inefficiencies in the 
tormer Soviet Union, the problem severity was much worse in Russia. Many research 


institutes faced extinction unless they could learn to compete in the new environment. 


Mr. John Kiser of International Scientific Products in San Jose, California, a 
consultant for the SDIO program, toured the former Soviet Union in 1988 to develop 
technology transfers to the U.S. During this tnp, he contacted the Kurchatov Institute, a 
developer of the TOPAZ-II system. (Malloy et al., 1994).This led to Russian participation in 
the Sixth Symposium on Space Nuclear Power Systems, held in Albuquerque, NM, in 
January 1989. At this symposium, the Russian delegation indicated their willingness to sell 
TOPAZ-type space power systems to the U.S. In a remarkable example of cooperation in 
the face of much institutional reluctance and bureaucratic inertia, the U.S. Government 
purchased two non-fueled TOPAZ-II systems and associated test and support hardware in 
1992. The SDIO funded investigation of the TOPAZ-II reactor system in support of the 
Nuclear Electric Propulsion Space Test Program (NEPSTP). This mission was intended to 
provided an early flight test with the inclusion of electric thrusters. SDIO was reconstituted 
as the Ballistic Missile Defense Organizaton (BMDO), resulting in NEPSTP program 
cancellation in 1994. The project continued as the Topaz International Program (TIP). 
(Wyant, 1994). Below is an image of a TOPAZ-II system at the NMERI facility (Figure 2). 





Figure 2. TOPAZ-II System at NMERI (From Malloy, 1994). 


The TIP is truly an international program. Scientists, engineers, and technicians 
from Russia, the United Kingdom, and France have conducted research at the Thermionic 
System Evaluation Test (ISET) facility in the NMERI complex in Albuquerque, NM 
(Fairchild et al. 1992). The New Mexico Alliance for Thermionics, which includes the 
University of New Mexico, Sandia National Laboratory, and Los Alamos National 
Laboratory, provides research support (Benke 1995). The BMDO Program provides 
funding through the, administration of the U.S. Air Force Phillips Laboratory. (The 
program may be shifted to the Defense Nuclear Agency.) 


The ‘TIP involves the testing and evaluation of TOPAZ-II system units, where 
nuclear fuel is replaced by tungsten heaters (Paramonov and Fl-Genk, 1994). The primary 
goal of this program is to gain understanding of the capabilities and hmitations of 


thermionic space nuclear power technology. Much of the basic thermionics research has 
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focused upon correlating insulator degradation to TFE lifetimes (Malloy et al., 1994). A 
stand alone TFE Test Rig allows TFE characteristic measurement. The test, monitoring, 
and control equipment of a Baikal Test Stand permitted performance evaluation of the 
TOPAZ-II power system. The TSET team has measured system parameters and 
characteristics in addition to various shock and vibration testing performed with reference 


to the MIL-STD-1540B flight qualification requirements. 





III. BASIC THEORY 


The TOPAZ-II system uses thermionic direct energy conversion with a nuclear 
fission reactor heat source. Direct energy (or static) conversion devices mainly use electrons 
as a working fluid rather than a vapor or gas characteristic of more conventional heat energy 
systems. This class of devices includes fuel cells, batteries, solar cells, thermoelectric and 
thermionic conversion mechanisms. A nuclear fission reactor is a device for producing heat 
from a controlled fission chain reaction. Fission is the process of splitting a heavy nucleus 


into smaller fragments and 1s accompanied by large energy release. 


A. THERMIONIC ENERGY CONVERSION 
1. General Characteristics 


A thermionic conversion system consists of two electrodes separated by a 
gap. The analogy of a heat engine with electrons a working fluid serves as a description of its 
characteristics. A “hot plate” or emitter 1s heated to a sufficiently high temperature that 
electrons are “boiled” off. These electrons travel to the significantly cooler collector where 
they “condense.” The gap between these electrodes is typically filled with an extremely low 


pressure gas or vapor. (Benke, 1995) 


2. Principles of Operation 


Electron emission from the emitter 1s similar to steam generation in heat 
engines. A thermionic converter may also be considered as a thermoelectric converter that 
uses plasma or vacuum as a conductor. The thermionic mechanism produces a potential 
difference between the two electrodes. When an external load is connected, this potential 
drives a current flow. By this process, part of the thermal energy that liberated the 


electrons is converted directly into electrical power. (Huang, 1988). 


3. Thermionic Efficiency 


The heat engine comparison is useful in understanding the factors which 


effect thermionic efficiency. The Camot efficiency, which 1s based upon the temperature 
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differences between source and sink, determines the maximum efficiency attamable. The 


main inefficiencies of the thermionic cycle are: 


a. Radiative Heat Transfer 


Heat transferred by radiation from the emitter to the collector 1s 
unavailable for energy conversion and is a loss. Heat rejection from the radiators is a loss 


inherent to the system. 

b. Thermal Losses 

Heat leakage through insulators and other undesirable paths reduces 
the efficiency of the cycle. 


© Back Emission 


Because the collector also has heat energy, some fraction of the 
collector electrons will have enough energy to deliver a current density to the emitter. This 
current reduces the overall forward electron flow. While the collector is at a lower 
temperature than the emitter, its lower work function results in the back emission losses 


described above. 


d. Space Charge Effect 


The emission of negatively charged electrons into the gap creates an 
overall space charge. A charge region creates an electrostatic field that inhibits the flow of 


additional electrons towards the collector. This effect limits the maximum current density. 


4. Thermionic Emission 


In metals and semiconductor materials under certain conditions, a layer of 
electrons is loosely bound in a region known as the Fermi cloud. The force binding an 


electron in this layer is given in Equation ( 1), below: 
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eR (1) 
l6z£x 


where. F = force (Newtons) 
e = electronic charge (Coulombs) 
E€ = permittivity of free space 


x = distance from surface (m) 


The work function, 9, describes the energy required to liberate the electron 


by overcoming the work function. The value of @ is unique to each material and has little 
temperature dependence in the normal range of thermionic conversion operation. Much 
thermionic material research concentrated on the development of desirable materials with 
good work functions and other acceptable properties in the high temperature environment 
within a converter’s enclosure. Current density from the emitter 1s a function of the rate at 
which electrons escape. The Richardson-Dushman equation, Equation ( 2 ), describes the 


relationship among the variables affecting emission current. 


ep 
Jo = RT exp) - 2) (25) 


where: J, = current density (amps / m?) 
R= material constant 
T — absolute temperature (K) 
e = electronic charge (coulombs) 
@ = work function (Volts) 


k = Boltzmann's constant 


This equation is useful in demonstrating the desirability of a high emitter 
temperature and the importance of carefully selecting materials (Angrist 1987). However, 
higher emitter temperatures also contribute to an increase of the radiation heat transfer as 


shown in the following equation, Equation ( 3): 
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Q- Aeol [uaa Di.) (3) 
where: O = radiative heat transfer (W) 

Á 7 area (m? ) 

€ — effective emittance 

o = Stefan - Boltzmann constant 

T — absolute temperature (K) 


This formula also stresses the importance of collector temperature. A lower 
collector temperature will increase radiation heat transfer from the emitter and reduces the 
amount of back emission. Higher collector temperatures are important in reducing the 
system radiator size, which also has a T° efficiency dependence at a cost however, of higher 
back emission losses. Optimum collector operating temperature requires a design trade-off 


study. (Angelo and Buden, 1985) 


5. Gas-filled Converters 


The strong impact of space charge build-up upon thermionic efficiency has 
direct consequences for converter design. Two methods of limiting this effect are 
minimizing gap size and introducing a plasma. The desired gap distance in a vacuum 
converter for useful current production is on the order of the Debye length. Debye length 
is the average distance traveled by an electron during the period of plasma isolation, and 1s 
on the order of 1 um (Rasor 1991). Because this represents a severely restrictive tolerance, 


vacuum converters have few practicable applications. 


A gas-filled (or plasma diode) thermionic converter uses space charge 
neutralization in order to increase gap width significantly (up to at most 0.5 mm) (Rasor, 
1991). This device continuously generates positive ions that neutralize the electrostatic field 


near the emitter by canceling the negative charge of the electrons. 


Cesium is the most commonly used gas used in plasma diodes because it has 
the lowest ionization potential of any element (Angelo and Buden, 1985). Cesium also 


promotes a low work function, @, which provides additional benefits. Cesium adsorbs to 


both the emitter and collector. Adsorption 1s the adhesion of a single layer of molecules to 
the surfaces of solid bodies or liquids with which they are in contact. This effect lowers the 
work functions of both electrodes, increasing overall converter efficiency. A liquid reservoir 
provides cesium through evaporation. Its production rate 1s dependent upon gas pressure 
and reservoir temperature. Efficient ion production requires emitter temperature 


approximately 3.6 times reservoir temperature (Angrist, 1987). 


6. Thermionic Modes of Operation 


The mode of operation in a thermionic converter 1s characterized by the 
mechanism of the neutralizing ion generation. In the unignited mode, surface ionization is 
the dominant production source. In the ignited or arc mode, inelastic collisions with 
electrons in the gap produce the majority of ions. The volume ionization in the ignited 
mode can generate a net positive space charge if the device is operated beyond the 
transition point. Collisions in the gap region cause ionization losses, creating a voltage 
deficiency 1”, Experimental data and theoretical investigation shows that minimum 1’, and 
maximum efficiency 1s achieved with small value of the product of cesium pressure (Pz, ) 
and electrode gap width (») (Angnst, 1987). Operationally, this performance sensitivity 


requires careful adjustment for converter optimization. 


B. FISSION REACTOR THEORY 
1. Fission 


Nuclear fission 1s the process in which the nucleus of an atom splits into two 
smaller nuclei. Fission 1s a spontaneous or induced event that usually requires massive atoms 
to occur. In a reactor, induced fission 1s primarily a result of the excitation of a nucleus with 
a neutron. Fission results in a large energy and radiation release, nuclear fission fragment 
formation, and the emission of more neutrons. A chain reaction occurs if these neutrons go 


on to induce additional fisstons. (Leachman, 1965). 


2° Neutron Energy 


The atomic absorption of an incident neutron forms a compound nucleus. 


Fissionable radioisotopes vary m the critical compound nucleus energy levels required to 
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induce fission. For example, ^?U is fissionable by slow neutrons whereas ^*U requires fast 
neutrons. The average energy of the neurons that maintain the chain reaction forms a basis 
for fission reactor classification. A thermal reactor operates with neutrons having a velocity 
distribution similar to that of materials at ambient temperatures. Intermediate or quasi- 
thermal reactors use higher neutrons energies. A fast reactor propagates the fission chain 


reaction with the high energy (~1 Mev) neutrons that are released from a fission event. 


Since a fission neutron 1s born with such high kinetic energy, it must be 
slowed to maintain a chain reaction a thermal or quasi-thermal reactor. A moderator 
matenal produces this slowing down by enabling scattering reactions. Fast reactors do not 


require moderation. 


S. Criticality 


Criticality is essentially a measure of the probability that neutrons produced 
in one generation will cause fissions in the next. A reactor 1s critical if, on the average, one 
neutron from a fission event induces a subsequent fission. In a critical reactor, both the 
fission rate and power are constant. The multiplication factor, &, describes the overall chain 
reaction condition. This factor 1s defined mathematically as the ratio of fissions produced in 


one generation to the number of fissions in the previous one (see Equation ( 4 )). 


k number of fissions in generation n + 1 
TT A TOS, AA (4) 


number of fissions in generation n 
In a critical reactor, k = 1. 1£ k < 1, a chain reaction will not be maintained 
and the reactor is said to be subcritical. If k > 1, the reactor 1s supercritical; power and 


fission rate increase. 


4. Neutron Life Cycle 


To maintain criticality, the chain reaction must produce one neutron for 
every neutron lost. Within each generation, a series of events defines the neutron life cycle. 
Although some fast fissions occur in all reactors, the majonty of the fissions in a thermal or 


quasi-thermal reactor require that a neutron be slowed without escaping. These neutrons are 
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absorbed within the reactor. An absorption into a non-fisstonable material (a nuclear 
poison), removes a neutron from the fission process. Additionally, only a certain percentage 
of the absorptions into fissionable material actually induce fission. These fissions provide 


the additional neutrons to renew the cycle into the next generation. 


Specific reactor. components are designed in order to change the 
probabilities that the average neutron will complete any step within the fission life cycle. A 
moderator increases the likelihood of neutron slowing down. Neutron absorpüve control 
devices (control drums or control rods) directly remove neutrons from this cycle. 
Additionally, many reactors are equipped with a reflector that minimizes leakage by 


deflecting neutrons back into the core where they may take place in the fission process. 


5. Reactor Control 


Reactor control is achieved through the manipulation of neutron absorbing 
devices within the reactor. These mechanisms use materials that have extremely high cross- 
sections for neutron absorptions. The most common control method govern criticality by 
changing the geometry of these absorbers with respect to neutron producing regions with 
the core. Since any factor that changes fission probabilities within the neutron Life Cycle 
will affect the chain reaction, other controlling mechanisms are possible. Routine regulation 
is required to start up and change reactor power and to overcome the effects of gradual fuel 
depletion and the build-up of poisonous fission product daughters. A reactor control 


system must also provide a means of emergency shutdown in order to ensure reactor safety. 


6. Delayed Neutrons 


Not all the neutrons produced dunng the fission. process are born 
immediately after fission. Some of the fission fragment daughters undergo further nuclear 
reactions that result in delayed neutron release. The half-life of delayed neutron production 
varies from about 1 sec to 1 minute. Since prompt neutrons are produced in about 10°" 
seconds after neutron absorpton, the effect of the delayed neutrons is to reduce the 
generation lifetime. This makes the process controllable by bringing the time scale into a 


range controllable by the reactor control system. 
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The fraction neutrons that are born delayed 1s the delayed neutron 
fraction, fJ. Since delayed neutrons are produced by decay processes, they are bom at a 
lower energy lever than prompt neutrons. The relatively slower delayed neutrons are less 
likely to leak out, and therefore more likely to be thermalized and cause fission. The 
effective delayed neutron fraction, B, accounts for this effect within reactor kinetics 


equations. 


7: Reactivity 


Reactivity, p, is the primary parameter used in reactor control theory. 
Reactivity is a measure of the departure from criticality. Reactivity is positive in a 
supercritical reactor, zero in a critical reactor, and negative when a reactor 1s subcritical. Any 
parameter change that increases criticality is said to add positive reactivity. Conversely, any 
event that tends to shut down the reactor adds negative reactivity. The primary reactivity 
control variables are control mechanism position and temperature feedback. Reactivity 1s 


mathematically defined with respect to criticality, as shown in Equation (5). 


EZ 


where: p = reactivity 
k = criticality 


Note: Reactivity is often reported in units of 10°, and is expressed 


with a '$' symbol. 


a. Control Drum Worth 


The TOPAZ-II reactor uses a control drum design in which sections 
of neutron absorbing material may be rotated. Turning the control drum moves this 
neutron absorber away from the center of the core. This decreases the probability that a 
neutron will reach the control drum, increasing the chance that fission will occur, adding 
positive reactivity. The reactivity of the control drums is with respect to position is the 


control drum worth. 
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b. Temperature Coefficient of Reactivity 


The temperature effect is especially important since it provides 
direct, immediate feedback to power changes. The parameter relating change in temperature 
to reactivity addition 1s the temperature coefficient of reactivity, a. A temperature change 
may add positive or negative reactivity, usually depending on the material involved. For the 
TOPAZ-II reactor, two effects are particularly important: Doppler broadening and Einstein 


oscillation. 


An increase in temperature 1s initially accompanied by an negative 
value of @, due to resonance absorption in **U. Resonance absorption is the process by 
which discrete neutron certain energy levels are preferentially captured as the neutrons slow 
down. These neutron captures do not result in fission, removing neutrons from the cycle. 
As the core temperature increases, a larger spectrum of energy levels are absorbed. This 


effect is known as Doppler broadening or spreading and adds negative reactivity. 


The Einstein oscillation effect is important in the Topaz-II 
moderator, which contains interstitial hydrogen. A higher moderator temperature increases 
the frequency of collisions and thermal interactions, which excites the hydrogen into higher 
energy levels. These higher energy levels have lower cross sections for capture and 
scattering, which results in more neutrons being available for the fission process. Therefore, 


the temperature increase results in the addition of positive reactivity, which contributes to a 


positive value of ar. 


8. Reactor Kinetics Equations 


The reactor kinetic equations relate power changes to the effects of delayed 
neutrons and reactivity changes. In this equation, the delayed neutrons are lumped in to N 
groups, based upon the decay constants of the precursors. These equations are valid for a 


small reactor with a uniform neutron flux distnbution. (Paramonov and El-Genk, 1994) 


2] 
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where: P = reactor power (W) 


p = effective delayed neutron fraction 
p(t) = reactivity 
A, = delayed neutron fraction of ith group 
£ = decay constant of the delayed neuttron precursor 
£” = prompt neutron lifetime 
C, = equivalent delayed neutron precursor concentration 


Q(t) 7 fission equivalent source power 


The initial conditions for Equations ( 6) and (7) are given below: 


P(6)= A (8) 
_ Ba 
C4) 


The reactivity input is calculated with Equation (9 ): 


E p.) PP Dar a 


external doppler temperature feedback 
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IV. TOPAZ II DESCRIPTION 


A. OVERVIEW 
1: Mission Requirements 


The TOPAZ-II space power system is a 6 kWe unit that was designed for 
geosynchonous missions with up to three year duration (Malloy et al, 1994). Figure 3 


shows an artist's conception of a TOPAZ-II spacecraft in a mission configuration. 





Figure 3. TOPAZ-II Mission Configuration. (From Schmidt et al., 1994). 


The general requirements for the TOPAZ-II design were: (Flight Safety 
Beam, 1992) 


1. System mass less than 1061 kg, excluding controller and power supply system. 
Electrical power 6 kWe, 27 volts, 3 year life span, 95% reliability. 


Shelf life of 10 years after manufacture. 


Toc er 9 


Reactor shall remain subcritical before reaching orbit. 


Reactor coolant must not freeze. 


(JI 


Among these requirements, the last two have directly influence the start-up 


control design. The interdiction against inadvertent criticality prior. achieving. mission 
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altitude 1s a nuclear satety concern. The prevention of coolant freezing limits the time 


available for the system to become self-sustaining. 


2, System Description 


The major TOPAZ-II subsystems are: the Reactor and combined 
thermionic converters, the Coolant System, Secondary systems, such as radiation shielding 


and cesium supply, the Power System, the Instrumentation and Control (I & C) System. 


The TOPAZ-II reactor provides electrical power with the generation of 
waste heat and radiation. Highly ennched uranium fuel heats the thermionic emitter, 
enabling electron flow. The liquid metal coolant system transfers waste heat to the radiator, 
limiting collector temperature. An electromagnetic (EM) pump provides the motive force 
for coolant flow. A radiation shield attached to the lower part of the reactor limits the 
neutron and gamma dose rate to the rest of the spacecraft. The cestum system supplies 
cesium (Cs) to the interelectrode gap, improving converter efficiency. A power supply 
system provides energy before the system is producing any usable electricity and regulates 
power production when the reactor 1s self-sustaining. The Instrumentation and Control 
(I & C) system monitors conditions, accomplishing start-up, operational control, and 
emergency shut down functions. 1 & C neutron and thermal detectors measure reactor 
power. The automatic control system (ACS) uses these signals to command control drive 
units which position control drums within the reactor subsystem. Figure 4 shows the 


locations of these subsystems. 
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Figure 4 TOPAZ-I System. (From Flight Safety Team, 1992). 





B. REACTOR SUBSYSTEM 
1. Overview 


The TOPAZ-II reactor is a small, zirconium hydride (ZrH, ss) moderated, 
epi-thermal design with a highly enriched *°U fuel. The Reactor Subsystem contains 37 
single-cell TFE’s that combine the fission heat source with the thermionic converters. The 
reactor subsystem also consists of a zirconium hydride (ZrH,,;) moderator, a beryllium 
reflector, 12 control drive mechanisms, and connections with the Cooling and Secondary 
Subsystems. Figure 5 shows a reactor top view, Figure 6 provides a side view. The TFE’s are 
cylindrical and fit into vertical slots within the moderator. Three of the TFE’s directly 
power the EM pump, the remaining 34 provide electricity to other TOPAZ-II and 
spacecraft loads. The TOPAZ-II design allows for fueling from the top of the reactor at the 


launch site, minimizing transportation safety concerns. (Schmidt et al., 1994) 
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Figure 5. Reactor Top View. (From Malloy et al., 1994). 
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Figure 6. Reactor Side View. 
2. Moderator 


The moderator material is ZrH,,, formed into block and enclosed in a 
pressurized stainless steel housing that also contains the working sections of the TFE’s, axial 
Be reflectors, and the connections to other systems. The axial reflectors are locate on the 
top and bottom of the moderator section. Interstitial hydrogen provides the high cross 
section for scattering required for fast neutron moderation. Therefore, maintenance of 
hydrogen pressure is essential for criticality. A pressurized mixture composed of 50% 
hydrogen and 50% carbon dioxide (CO,) within the moderator inhibits hydrogen release 


and improves heat transfer. The reactor housekceping telemetry reports this moderator gas 
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pressure. Reactor materials and coatings are specifically designed to minimize hydrogen 


diffusion. (Voss, 1994) 


3. Thermionic Fuel Elements 


Fach TFE has an inner fuel section enclosed by the emitter, gap and 
collector of the thermionic converter. An insulating layer separates the collector from a 


stainless steel outer jacket. Figure 7 shows the TFE with its associated spacers and 






insulators. 
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Figure 7. Thermionic Fuel Element (After Voss, 1994). 


a. Fuel Section 


The fuel is composed of 96% *°U enriched uranium oxide (UO,) 
formed into cylindrical fuel pellets with a hollow center. This center section is the release 
path for radioactive fission products. Beryllium (Be) and Berylltum oxide (BeO) reflector 
pellets are positioned on top and bottom of the fuel stack. A heater assembly replaces the 
fuel section for non-nuclear testing. Nominal fuel temperature is between —1500 C to 
1650 C. (Schmidt et al., 1994). 


b. Thermionic Section 


The emitter is a monocrystal of molybdenum with 3% niobium 
added for strength. A layer of tungsten formed by chemical vapor deposition establishes the 
desired emitter work function. This tungsten layer is in the form of the "^W isotope, 
minimizing parasitic neutron capture. Helium (He) fills the 0.5 mm width interelectrode gap 


prior to start-up, providing superior heat transfer characteristics. Cesium replaces He after 
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E 


the Cs supply system is placed in operation. The collector is formed from polycrystalline 
molybdenum with a plasma sprayed alumina (AI,0;) insulator outer coating. The thermionic 
section transforms approximately twenty percent of the heat energy from the emitter into 
electron flow. The collector temperature is approximately 627 C (900 K) during operation. 
A helium gap, located between the collector insulator and the inner coolant tube, which 
improves heat transfer while maintaining electrical insulation. Figure 8 provides a cross 


section of a TFE working section. (Voss, 1994) 
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Figure 8. TFE Working Section. (After Flight Safety Team, 1992). 
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4. Radial Reflector and Control Drums 


A beryllium (Be) reflector encloses the moderator section. The Be material is 
divided between static Be inserts and the twelve rotating drum assemblies (See Figure 5). 
Two stainless steel bands hold the radial reflectors in place against spring pressure. An 
electrical signal from the automatic control system can initiate reflector release, providing 
emergency shut-down ability. High temperatures or mechanical shock equivalent to a 50 m 
fall will also split the bands. The Be control drums have a section of boron carbide (B,C) 
that spans 120 degrees. Boron has a very high cross section for neutron absorption. A 
control drive system positions nine control drums nearly simultaneously. An independent 


mechanism positions three of the control drums, the safety drums, during start-up. 


E. HEAT REJECTION SYSTEM 
1. Overview 


The sodium potassium (NaK) coolant system removes waste heat from the 
reactor. The schematic diagram provided below (Figure 9), shows the connections between 
coolant system components: the direct current powered EM pump, a volume accumulator, 
a radiator, electrical heaters, piping, pressure and temperature monitoring equipment. The 
EM pump provides the motive force for coolant flow. The volume accumulator allows for 
thermal expansion and contraction. The electric heaters are used during the pre-launch 
preparation phase to prevent NaK freezing. Waste energy from the thermionic process 
heats coolant from 743 K to 843 K dunng its single pass through the core to the upper 
plenum. Coolant then flows to the radiator, where heat is rejected to space. The EM pump 
returns the coolant to the lower inlet plenum, where the cycle repeats. (Paramonov and 


El-Genk, 1994) 
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Figure 9. Primary Coolant System. 


The design requirement that coolant freezing be prevented is critical to 
system startup operation. Since NaK freezes at -11 C, the design minimum temperature is 
-5C. Prior to launch, electric heaters maintain coolant temperature. Prior to startup, 
intermittent the operation of the EM circulates coolant, equalizing NaK temperatures using 
residual heat. A thermal blanket covers the radiator in order to limit heat losses. This 
blanket is ejected shortly after the reactor reaches the point of adding heat. The cover was a 


requirement for the geosynchonous orbit of the Soviet flight test. A U.S. mission with a 
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lower operational orbit may not need a blanket due to the reduced time from launch to 


reactor startup. (Flight Safety Team, 1992) 


2: Electromagnetic Pump 


The EM pump is a DC powered unit that operates on the Faraday principle. 
Since the Nak coolant 1s a liquid metal, it 1s a good electrical conductor. An electromotive 
force is generated perpendicular to an electric current and a magnetic field. During normal 
operating conditions, three TFE’s connected in parallel provide EM pump power. Startup 
batteries power the EM pump intermittently prior to reactor startup. During startup, battery 
current limits the flow rate to 25% (~0.4 kg/s). When coolant temperature in the lower 
collector reaches 235 C, the battery current switches to maximum. Flow rate increases from 
50% (~0.8 kg/s) to 100% (~1.3 kg/s) as reactor power increases and the TFE’s provide 
useful voltage. (Kwok, 1993) 


D. SECONDARY SYSTEMS 
1. Radiation Shielding 


A lithium hydride, stainless steel encased, radiation shield attenuates neutron 
and gamma radiation. The shield 1s located behind the reactor in order to minimize dose 
rate to the payload and other TOPAZ-II components. Penetrations in this shield are angled, 
which reduces radiation streaming. Since radiation levels outside the reactor are significantly 
above background, an eight meter boom will separate the TOPAZ-II SNPS from the 
payload and spacecraft bus, as previously shown in Figure 3 (Flight Safety Team, 1992). At 


this distance, radiation spreading sufficiently minimizes the effect upon the payload. 


2 Cesium Supply System 


The Cesium Supply System (CSS) provides Cs to the interelectrode gap in 
order to increase thermionic efficiency. Figure 10 shows the overall system diagram. Cs 
condenses on a small radiator within the Cs regulator (see Figure 11). A “wick” formed of 
wound stainless steel wires transports Cs from the radiator. Cs vaporizes as it moves along 
the wick, establishing a relatively constant supply during normal operating conditions. A 


throttle valve establishes a nominal Cs pressure of 2 Torr to the reactor. A release line 
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provides a bleed-off path through which —0.5 g/day of Cs and impurities vent into space. 
This process maintains the required Cs purity of 99.9999, (Flight Safety Team, 1992). Since 
only —1 kg of Cs is supplied, the loss of Cs ts a limiting factor for mission duration (Schmidt 
et al., 1994). 
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Figure 10. Cesium Supply System. 


Helium (He) gas fills the interelectrode gap prior to startup. He is used 
because of its non-reactivity and excellent thermal conductivity. As the CSS heats up, Cs 
mixes with He. During the startup sequence, a puncture valve (V3), opens the release path 


into space. Cs then replaces He, which is preferentially vented due to its high diffusion rate. 
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Figure 11. Cesium Regulator and Reservoir (After Morris et al., 1994). 


E. INSTRUMENTATION AND CONTROL 

The Instrumentation and Control (I & C) system provides for the remote monitoring, 
control, and reporting of system status. The I & C subcomponents include the measuring 
instruments and the Automatic Control System (ACS). 

1. Instrumentation 

TOPAZ-II measuring instruments include devices for monitoring neutron 
power, temperatures, pressures, and control drum positions. 

a. Neutron Flux Measurement 


A measurement of neutron flux provides indication of reactor 
power. Two ionization chambers (IC’s) containing uranium coated plates are located 


adjacent to the reactor. Nuclear reactions of the neutrons with the uranium produce 
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ionization, which. enables current flow. When neutron flux 1s at a level corresponding to 


100% power, the IC output 1s approximately 1 ma. 


The ACS uses the highest of the two IC currents as its control input 
for reactor power level. Since the Russian designed IC’s are insensitive to power levels 
below ~6% power, the U.S. design will include two to four hafnium or gadolinium Self- 
Powered Neutron Detectors (SPND's). Not only do SPND's provide accurate indications 
of low neutron power, they can be used as a reliable alternate indication at higher power 
levels. The SPND signal would be suitable for use as a safety monitoring system input. 
(Flight Safety Team, 1992) 


b. Temperature Measurement 


The 18 C system uses two types of temperature monitoring: 
Chromel-Alumel thermocouples (T/C’s) and platinum resistance temperature detectors. 
Paired T/C’s measure reactor coolant inlet and outlet temperatures inside the pipe wall. 
These T/C’s have a six second response time, which may limit their utility as an input to a 
reactor control scheme. The RTD's measure coolant temperatures at the radiator, control 


EM pump operation during startup, and allow for system monitoring prior to launch. 


ci Pressure Measurement 


Figure 12 presents a diagram for the pressure sensors in the I & C 


system. Pressure sensors monitor the following parameters: 


1. Moderator gas pressure. 

TFE gap gas pressure. 

Shield gap gas pressure. 

NaK volume accumulator gas pressure. 


Startup battery electrolytic pressure. 
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Control drum assembly pressure. 
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Figure 12. Pressure Sensor. (Flight Safety Team, 1992). 


The pressure signals are used for telemetry purposes only, they do 
not provide any control function inputs. Since the volume accumulator pressure would 


provide indication of a coolant leak, it may be useful as a safety system input. 


d. Control Drum Position Determination 


There are three means of control drum position monitoring. The 
primary method uses the control drum position sensor. The sensor provides an analog 
measurement of position indication. A roller assembly within the control drum unit 
positions a rotor element. This rotor varies the coupling of a 1 kHz signal in a transducer, 
causing an amplitude change in signal coil output. The backup method counts control 
system signal pulses to the drive motor, which provides a more accurate measure of drum 
position as long as the control drum motor 1s operating normally. Additionally, limit 


switches indicate drum travel has reached the upper (154°) or lower (0°) bound of travel. 
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2: Control System 
a. Overview 


The ACS components are the Telemetry and Command System 
(TCS), the Automatic Regulation System (ARS), and the Electnc Power Supply System 
(PSS). The TCS processes commands, relays telemetry, and sequences startup actions. The 
ARS receives inputs from the instrumentation and generates a control signal to position the 
control drum assembly. It also conducts shutdown and protective control actions. ARS 
design requirements specify a thermal control range of 0 to 150% power with an accuracy of 
+700 W. The PSS regulates power supplied to the spacecraft bus, maintaining bus voltage at 
27+0.8 volts. John Hopkins Applied Physics Laboratory (APL) has designed a modified PSS 
with a lower weight and without the gas pressurized coolant requirements of the Russian 


system. 


b. Control Drum Drive System 


Each control drum contains a B,C neutron absorbing poison plate. 
The positioning of each poison plate with respect to the main core controls reactor 
criticality. The control drum drive system rotates the nine control drums. The maximum 
angle range 1s from 0° to 180° at speeds of 1.4°/s. The drive assembly is sealed in an argon 
gas pressurized enclosure since the mechanism 1s not designed to operate in a vacuum. The 
assembly consists of a stepper motor, drive train, main drum, position indication devices, 


temperature sensors, and pressure sensors. 


A signal from the ARS powers the stepper motor. A gas tght 
bellows transfers torque from the reduction gear to the drive train. The drive shaft backlash 
is 1° to the main drum. This drum has slots through which cams connected to the bottom 


of the control drums are linked. (Flight Safety Team, 1992) 


The control drum cam and main drum slot design enables a very fine 
control of thermal power. Each control drum will only rotate once the slot edge comes in 


contact with the cam. On a drum reversal, only one drum will move initially, followed by 
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other drums as more slot edges make contact. This feature establishes a planned hysteresis 


of 9° in moving the control drums. (Voss, 1994) 


3. Protection System 


The protection subsystems within the TOPAZ-II design are the Safety 


Drum Drive System and the Reflector Ejection system. 


a. Safety Drum System 


Three drums with a combined reactivity worth of -$2.0 are held in 
the “IN” position prior to reactor startup. This provides sufficient shutdown margin to 
prevent an inadvertent criticality accident in most probable accident scenarios not involving 
water flooding or immersion. The NEP Space Test Program Preliminary Nuclear Safety Assessment 
recommends addition safeguards in order to meet U.S. design criteria for these additional 
possibilities. Each safety drum is connected to a direct current drive motor that is controlled 
by the TCS. After the TCS initiates safety drum rotation dunng startup, the drums rotate 
until a limit switch cuts off power at ~180°. (Flight Safety Team, 1992) 


b. Reflector Ejection System 


The reflector system provides a quick means of initiating emergency 
reactor shutdown. The reflector segments and drums are held against spring pressure by a 
stainless steel band with fusible links. A high current signal from the TCS will melt the links 
in ~0.25 seconds, releasing the bands, which causes reflector disassembly. The loss of the 


reflector reactivity will leave the reactor subcritical. 


The Russian safety system design included three pairs of events that 


would initiate reflector ejection: 


1. Receipt of 2 ejection commands from the ground, in senes. 
2. Receipt of 2 emergency rocket shut off signals. 


3. Two launch vehicle emergency signals. 
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The Russian design did not include any provision for reflector 
ejection once the system was operating in orbit. U.S. design philosophy would require this 
protection, especially in view of the much lower altitude mission profiles envisioned for a 


TOPAZ-I flight test. 
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Y. SAFETY AND FUNCTIONAL REQUIREMENTS 


A. NEPSTEP MISSION 
1. Overview 


The Nuclear Electric Propulsion Space Test Program (NEPSTP) mission 
was an outgrowth of Amenca’s Space Exploration Initiative, developed dunng the Bush 
administration (Synthesis Group, 1991). Since current policy precludes any U.S. reactor 
launches, NEPSTP development remains suspended, and exact mission details are 
undetermined. The onginal flight evaluation mission plan called for a launch into 1000 km 
to 2000 km transfer orbit, where the TOPAZ-II reactor system startup would begin. Due to 
concems for the system's impact on the space environment, this transit orbit altitude was 
raised to 5250 km, minimizing interference with the Cosmic Ray Observer mission. After a 
series of tests, electrical thrusters using xenon propellant will increase the altitude until the 


spacecraft attains a 36000 km circular high Earth orbit (HEO). 


2 Purpose 


The primary mission goal is to evaluate the use of nuclear electric propulsion 
for orbital transfer applications, with particular emphasis placed upon measurmg the self- 
induced environment produced by the system when operating. The secondary goal 1s to test 


diverse electric thrusters and evaluate the performance of the TOPAZ-II system. 


3: Mission Details 


The mission 1s divided into three main phases: pre-operational, operational, 
and post-operational. Figure 13 shows a diagram of these phases. The pre-operational phase 
includes launch and checkout. During the operational phase, nuclear electric propulsion 
raises the spacecraft to HEO. Reactor shutdown and system venting occurs in the post- 


operational phase. 
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Figure 13. NEP Mission Phases & Scenarios. (From Marshal, 1994). 
a. Pre-operational Phase 


The planned NEPSTP mission launch site is the Cape Canaveral Air 
Force Station, in Cape Canaveral, Florida.. The spacecraft mass will be approximately 
3500 kg, which will require a medium sized launch vehicle to reach transfer orbit. This initial 
circular orbit will have an inclination angle of 28.5° at a minimum altitude of 5250 km. Both 
ground and launch vehicle signals will be capable of initiating emergency reflector 


disassembly in order to minimize the probability of inadvertent criticality upon vehicle 
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reentry. Prior to launch, final assembly, nuclear fueling, and systems checkout will venfy 


readiness. 


b. Operational 


Once the spacecraft has been placed into transfer orbit, the payload 
boom will extend from the TOP AZ-II reactor system. Since the Nak coolant is kept from 
freezing by residual heat, reactor startup must be initiated within. about eight hours of 
launch. Within one hour from extension, the reactor should be self-sustaining, with the 
TFE system powenmng all loads and performing battery charging. Following a series of tests, 
the nuclear electric propulsion system will be activated. The thrust will be directed along the 
vehicle's velocity vector, increasing alttude while maintaining an. essentially circular. orbit. 
The TOPAZ-II system will take an estimated 500 days to reach about 3600 km altitude. 
(Flight Safety Team, 1992) 


For most of the operational period, the spacecraft will be unable to 
communicate with the mission control center. This has important implications for control 
system design since the spacecraft must be capable of autonomous control. This requires 
the addition of supervisory control and a means for an emergency shutdown, features not 
found in the Russian design. Additionally, provisions must be made for automatic shutdown 


if contact with the ground station 1s lost for an extended period of time. 


c. Post-operational 


In accordance with requirements to minimize radiological 
consequences of the NEPSTP, safe disposal requirements are an essential part of mission 
planning. Once the mission 1s complete, the reactor must be shut down with systems 
vented. Disposal requires that the system be in a "sufficiently high” orbital alutude. A 
sufficiently high orbit is one in which the expected orbit lifetime is long enough to ensure 
that the curie content of the core decays to the intrinsic actinide levels of the fuel prior to 


bum-up. 
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B. SAFETY 
1. Philosophy 


U.S. and Russian design philosophies differ. The Russian design was tailored 
for a specific mission involving initial criticality at geosynchonous orbit. The U.S. Topaz-II 
Flight Safety Team performed additional preliminary safety analysis since the U.S. mission 
scenario includes possible accidents not viable within the Russian flight profile. Figure 14 
shows profiles of some pre-operational accidents. Russian design emphasized operational 
reliability over the risk of premature shutdown. No provisions were made for automatic 
shutdown at the end of the mission. The reactor would eventually shutdown due to the loss 
of reactivity caused by hydrogen leakage and other expected phenomena. Furthermore, 
reflector ejection, the primary rapid means of reactor shutdown, 1s inoperable in the Russian 
design after the spacecraft is released from the launch vehicle. The only operational 
safeguard is the temperature regulator, which limits coolant temperature to less than 600 C. 
(Flight Safety Team, 1992) 
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Figure 14. Pre-operational Accident Scenarios. (From Marshal, 1994). 
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2. Safety Requirements 


The primary goal of the Topaz-II satety policy 1s to reduce risks to as low as 
reasonably achievable. This policy has two main concerns; safety and environmental 
protection. Safety requirements are a fundamental consideration for any nuclear program. 
Many safety provisions are legislated in U.S. and international law. Since a reactor accident 
involves the potential release of radioactive materials, a space nuclear power mission must 
comply with U.S. environmental law and obligations under Article IX of the Outer-Space 
Treaty, which provides for international protection of the outer-space environment (Larson 
et al, 1992). Implementation of these requirements has provisions that fall into the following 


categories: 


a. Reactor Startup 


The reactor shall not be operated at significant power levels until the 
spacecraft reaches a sufficiently high orbit. The basis for this provision. that the curie 
content of the core is based upon the time-history profile of the activation processes of the 
fission reaction. Pre-orbital operation of the reactor could expose the ground crew and any 


emergency impact site to undue risk. 


b. Inadvertent Criticality 


A reactor relies upon a specific geometry and placement of materials 
to maintain a sustainable chain reaction. All plausible accident scenarios must be analyzed to 
ensure that inadvertent criticality does not result. Impact induced criticality would be 
uncontrollable since control and monitoring systems would also fail. A subcritical reactor 


will not increase its radioactive inventory and emits far less gamma and neutron radiation. 


Initial results of Monte-Carlo-based neutronic analysis show that the 
reactor would be supercritical following water flooding and immersion with intact radial 
reflectors. A wet sand scenario produced the worst-case configuration. Therefore, in order 
to prevent inadvertent criticality, the U.S. implementation of the design requires 


modifications. The current proposal plans for keeping four fuel elements external to the 


core until the initiation of reactor startup. This change will ensure that an accidental critical 


configuration would not be attained. 


3. Radiological Release 


The release of radioactive material during operations and following any 
probable accident scenarios must be minimized. This provision reduces any 1mpact upon 
the Earth environment and decreases the potential interference with other space missions. 
Consideration of the Cosmic Ray Observer mission led to an increase in the operational 


orbit altitude. 


4. Disposal 


The TOPAZ-II design and mission. requirements stipulates provisions for 
safe disposal. System operation shall be limited to a sufficiently high altitude for on-orbit 
disposal. Unlike the initial Russian design, the control system will assure reactor shutdown at 


the end of mission life. (Flight Safety Team, 1992) 


5, Re-entry 


Prior to initial high power operation, the TOPAZ-II system does not 
present a significant radiological hazard. Safe re-entry considerations rely upon following 
requirement for reactor startup and the prevention of inadvertent criticality. The system 


mechanical design provides for minimal debris dispersal upon impact. 
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VI. SYSTEM STARTUP 


A. MISSION PHASES 


The following sequence 1s based on the NEPSTP mission profile with the incorporation of 


proposed U.S. design changes. 


1. Pre-Launch 


Spacecraft assembly includes reactor fueling and low power testing. Prior to 
launch, an umbilical monitors system operation, charges the startup battery, and provides 
power for Nak coolant heating. A T,,,, regulator maimtains the liquid metal coolant at about 
100 C. The minimum temperature for launch is 40 C, which is based upon maintaining 
enough heat to provide an eight hour margin to prevent NaK freezing. The Automatic 


Control System (ACS) is also energized. 


2. Launch 


Immediately after launch, the reflector eyector system is enabled. A series of 
two signals (an “allow” signal and an “eject” signal) sent from the launch vehicle or ground 
control can start the reflector ejecnon sequence. The reflector ejection initiation signal ts a 
0.1 amp pulse of 0.1 second duration that closes an eject contact. Closure of this contact 
sends high current from the startup battery through the ejection band locks, which melt 
within 0.25 seconds, allowing spring force to cause reflector ejection. (Flight Safety Team, 


1992) 


During the orbital transfer phase, NaK coolant temperature ts maintained by 
periodic cycling of the electromagnetic (EM) pump. The pump circulates the residual heat 
energy, leveling NaK temperatures. Estimated battery capacity allows for approximately 
eleven cycles of 90 second pump operation prior to TFE power production. Figure 15 gives 
a graphical representation of the results of a simulation of temperature regulation action. 


This analysis agrees with the estimate of eleven cycles. Note that most of the pump 
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operation occurs during the final hours of this phase. Shortening the orbital transfer time 


may significantly reduce the required startup battery capacity. 


NaK Post-Launch Temperature Profile 
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Figure 15. Post Launch Tmin Regulation. 
3. On-orbit 


After the spacecraft has been placed in a stable minimum altitude orbit, the 
boom deploys, separating the TOPAZ-II from the payload and control system. A series of 


tests will ensure that the system is ready for startup. 


4. Startup Phase 


The Russian design would perform an automatic reactor start following 
boom deployment, with a ground command signal as a backup. John Hopkins Applied 
Physics Laboratory (APL) recommends enabling startup with ground command only. This 
modification brings the system into closer compliance with U.S. safety. The two startup 


signals are “Unlock Startup” and “Startup.” 
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B. 


STARTUP SEQUENCE 


Startup Limits 


Thermal limits form the basis for the operational limits during startup. The 


peak core temperature (Ta) is limited to a maximum of 2000 C. The temperature limits for 


the reactor moderator and vessel are 800 C. Table 1 presents the resulting operational 


restrictions. Operation outside these limits will result in undesired changes in material 


properties and place excessive mechanical stress upon the system. 


wise linear control law. The design criterion for the controller are to maintain thermal 












Operating Power 
Range (kW 


115 Steady State 600 C 
Operation 


Table 1. Reactor Power Startup Limits. 





Figure 16 shows the time history of the control band formed by this piece- 


power between 0 to 150% of maximum nominal thermal power of 115 kW. 
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Figure 16. Reactor Power Startup Control Band 
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2: Startup Sequence 


Receipt of the “Unlock Startup” signal tums on the Automatic Regulating 
System (ARS). This enables the neutron power and coolant temperature limit 
regulation control of the control drum system in accordance with the startup 
limits. Thermal cover ejection 1s also enabled. 

The four withdrawn fuel elements are inserted. Since the mechanism that would 
restrain and insert these elements is not yet designed, the time required for this 
step is uncertain. However, it is probably on the order of ten seconds. 

The three Safety drums are rotated outwards, providing +$2.06k worth of 
reactivity from the shutdown value of -$6.0 5k (Paramonov and El-Genk, 1994). 
This step requires approximately 8.5 seconds. The current safety drum drive 
motor is not wired for reverse rotation, which would provide a redundant 
backup safety feature. Since these motors are not designed to withstand the 
harsh radiation environment of an operating reactor, rotation inward could only 
be possible prior to operation at power.. Therefore, the usefulness of wiring the 
motor for reverse operation requires analysis of probable startup casualties. The 
capacity for reverse safety drum operation after full reactor startup requires 
complete redesign of the motor. 

Control drums positioned to 154° and back to 145° in accordance with preset 
limits. Although this step uses open loop control, the ARS regulation is 
operational and would initiate control drum inward rotation if a set-point is 
reached. At this stage, ARS action could only occur following some casualty. 
This regulation is the primary reactor safety response for analysis purposes. The 
initial control drum movement will provide excess reactivity $0.30 5k at 0 C to 
$0.70 5k at 70 C. The reactor is now supercritical, with power increasing with 
period of 20 to 4 seconds, respectively. 

At Five percent (5%) reactor power, the Minimal Control Level (MCL) for the 
Russian designed neutron power detectors 1s reached. Minimal Control Level 1s 


that power level in which the neutron power level ts sufficiently above 
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11. 


17 


13. 


background noise to provide a reliable reading. Inclusion of Self-Powered 
Neutron Detectors (SPND’s) will provide a neutron signal over a wide response 
range from a low level neutron reading up to approximately eight percent (8%) 
reactor power. The SPND's wil also provide useful telemetry for studying 
TOPAZ-II startup operation. 

At one to three percent above MCL, approximately 6-9% power, the neutron 
power regulator switches on, shifting the ARS to closed loop control. Primary 
feedback control is provided by neutron power rate regulation; the coolant 
temperature limit system provides a safety backup. 

The MCL signal also turns on the electromagnetic (EM) pump, establishing 
sodium-potassium (Nak) coolant flow at 25% of full flow. 

When coolant temperature reaches 100 C, the thermal cover ejected. 

After the coolant has heated to 230 C, the control system provides maximum 
battery power the EM pump, increasing flow to approximately 50%. The EM 
pump 1s supplied in parallel by the load current produced by its three dedicated 


thermionic fuel elements (ITFE’s). This current is negligible at this point. 


. At a nominal neutron power level of 110 kW, the upper limit, the neutron 


regulator stabilizes power. 

Useful thermionic power is generated as cesium (Cs) pressure increases. The Cs 
system is warmed up with the NaK system by way of a heat exchanger. The 
power from the TFE's is converted into heat by a series of load resistors. 

Once a TFE load current of 60 amps is reached, the cesium vent valve opens. 
The Helium-Cesium (He-Cs) mixture in the TFE interelectrode gap is replaced 
by Cs as the He vents out. 

The power supply system (PSS) regulator switches on when output power 
reaches a current level 2 100 amps and a voltage level of > 27 volts. The PSS 
establishes the neutron power set-point, within design limits. The PSS directs 


TFE power to the payload. 


. Once the PSS is operating, the startup battery switches are opened and the 


battery electrolyte is vented into space. Electromagnetic pump power 1s supplied 


solely by its three TFE's. 


S 


15. The startup sequence 1s complete once the reactor system 1s producing steady, 
reliable power. The Russian mission profile allowed for 70 minutes for startup. 
The time required by the Nuclear Electric Propulsion Space Test Program 
(NEPSTP) mission will be different, since 1t will be commenced from a different 
altitude, and only when the spacecraft is in communications with the ground 
station. 


E SYSTEMS TESTING 


In addition to the non-nuclear testing within the Topaz International Program, TIP, 
the system has been extensively analyzed within the Russian test program. Among the test 
results that have a direct effect upon startup control are: the Russian Plant 82 Nuclear 
Ground Test, in which a loss of coolant accident (LOCA) occurred, and the cesium 
pressure oscillations observed at testing at the New Mexico Engineering Research Institute 


(NMERI). 


1. Loss of Coolant Accident 


A loss of coolant accident occurred after approximately one year of nuclear 
testing at the Plant 82 ground site. A NaK leak developed from a crack in a weld on an EM 
pump discharge line (Flight Safety Team, 1992). The initial symptoms of this event were an 
in increase in the temperature rise across the core of approximately 15 C and current and 
voltage oscillanons. The coolant temperature and pressure conditions led to boiling within 
the coolant channels, causing power oscillations of about 4 kW,,. The loss of flow casualty 
led to a loss of all coolant as the NaK was released into the vacuum test chamber. The 
resulting increases in core temperatures led to the deformation of structural materials and 
the development of cracks within the ZrH,,, reflectors. The release of hydrogen following 
this damage added enough negative reactivity to shutdown the reactor, even though the 
operators had attempted to maintain criticality by manually rotating all control drums 


outward. 


This LOCA has several implications for the successful implementation of a 
ACS design. The system demonstrated the robustness of its design. This accident did not 
develop into a radiological hazard even though operator action in attempting to maintain 


power was only compounding the casualty. The TOPAZ-II has a demonstrated structural 
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safety margin with adequate provisions for heat transfer without coolant to prevent serious 
damage to the TFE's or the reactor vessel (Flight Safety Team, 1992). Both the U.S. 
designed ACS and the operating procedures must preclude continued operation following a 
ROCA. 


This accident also highlighted limitations of the temperature sensor package 
planned for the flight system. There are no in core temperature sensors. The ground system 
temperature sensors are located on the inlet and outlet coolant pipes. The flight system has 
sensors on the outlet pipes only. Due to the leak location, these sensors did not detect the 
coolant temperature increase following the LOCA. Therefore, the temperature limit 
regulation cannot be relied upon as the sole indication of this type of casualty. Without 
further supervisory safety monitoring, the ARS would react to the loss of current produced 
by a LOCA by attempting to increase reactor power. The lessons learned from this incident 


rust be analyzed and incorporated into control system design. 


2: Cesium Oscillations 


The first TOPAZ-II unit tested as part of the TSET program , the V-71 
prototype, experienced periodic oscillations in its output current. The second unit tested, 
Ya-21U, also exhibited wide swings in power during startup These oscillatons were 
enigmatic since the Russian tests did not reveal this behavior (Morris, et al., 1994). The test 
facility at NMERI includes a Russian supplied “Baikal” vacuum test stand with U.S. control 
and monitoring modificanons (Wold, 1994). The root cause for these oscillanons has been 
traced to instabilities in cesium pressure (Moris, et al., 1994). The cesium system relies upon 
a careful balance between pressures and temperatures to supply a constant Cs flow. The 
presence of cold spots in the Cs system, particularly during startup, disrupts this equilibrium. 
This effect 1s significant from a controls system standpoint since the tracking and stabilizing 


these oscillations is problematc. 


The first TOPAZ-II unit tested as part of the TSET program , the V-71 
prototype, experienced periodic oscillations in its output current. The second unit tested, 
Ya-21U, also exhibited wide swings in power during startup These oscillatons were 


enigmatic since the Russian tests did not reveal this behavior (Morris, et al., 1994). The test 
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facility at NMERI includes a Russian supplied “Baikal” vacuum test stand with U.S. control 
and monitoring modifications (Wold, 1994). The root cause for these oscillations has been 
traced to instabilities in cesium pressure (Morris, et al., 1994). The cesium system relies upon 
a careful balance between pressures and temperatures to supply a constant Cs flow. The 
presence of cold spots in the Cs system, particularly during startup, disrupts this equilibrium. 
This effect 1s significant from a controls system standpomt since the tracking and stabilizing 


these oscillations ts problematic. 
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VII. 


STARTUP CONTROL 


A. AUTOMATIC REGULATION SYSTEM 


Automatic Regulation System (ARS) functions are divided between three separate 
controllers, as shown in Figure 17. The following parameters apply to the Russian ARS 
design The Supply Current Regulator establishes the reactor power setpoint based upon 
TFE supply current, with a £3.5 amp deadband. The Reactor Power Regulator compares 
reactor power as measured by the neutron detectors to the this setpoint and sends 
commands to the Drum Drive Controller. The Reactor Power Regulator deadband is £250 
watt. The Drum Drive Controller directs operation of the control drum motor. The 
purpose of the coolant temperature limiter is supervisory. If reactor coolant outlet 
temperature exceeds 585 C, the neutron power setpoint cannot be increased. If the coolant 
temperature reaches 600 C, the temperature limiter lowers the setpoint, resulting in inward 


motion of the control drums. The Telemetry and Command System (TCS) provides overall 


system control. 
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Figure 17. Control System Block Diagram. 
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This system has two modes of regulation: Startup Mode and Operating Mode. In 
startup mode, the reactor power follows the preset control band shown in Figure 16 
(above). After startup, the regulator shifts to Operating Mode and the Supply Current 
Regulator establishes the ordered neutron power level. The coolant temperature limiter 


provides safety override. 


B. SYSTEM MODELING 
1. Current Models 


There are many models available for simulating the TOPAZ-II system. 
Three models that are of particular interest for startup control design are the Thermionic 
Transient Analysis Model (ITI AM), the Reactor System Dynamic Simulator, and the Sliding 
Mode Control Model. Figure 18 shows the common elements within these simulations. The 
neutronic model includes the reactor kinetics equations. The thermal model accounts for 
the temperature changes throughout the system and provides reactivity feedback to the 
neutronic model. The electrical model includes the thermionic conversion processes, which 
affect system temperature through electrical energy transfer. The control system model 


includes the control law and the control drum drive response. 
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Figure 18. TOPAZ-I Modeling. 
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a. Comparison of Simulations 


Analysis of these three mathematical models 1s based upon 
comparison of their component models and an examination of the differences among their 
results. This permits development of a hybrid control model that incorporates desired 


features from each simulation. 


D TITAM 


The Thermionic Transient Analysis Model is programmed with 
modular components that sımulate the neutronic, thermionic, thermohydraulic, and 
electrical behavior of the TOPAZ-II system. This model was developed as part of the 
Topaz International Project research. The program simulates transient and steady-state 
operation with either fission or electric heating. The equations within this model are based 
upon a complete, mgorous application of the physics principles behind the TOPAZ-II 
technology. The TITAM results show excellent agreement with both the Russian 
simulations and measurements taken at the TSET facility. (Paramonov and El-Genk, 1994). 
Since a benchmark of this simulator has been successfully performed, the hybrid control 
model will use its equations unless there 1s sufficient justification for another choice. The 
authors also generously provided a modified version of the TITAM simulation that would 
allow analysis of step power transients and the determination of equilibrium conditions at 


different power levels. 


c. Reactor System Real-time Dynamic Simulator 


The TOPAZ-II Reactor System Real-time Dynamic Simulator is an 
IBM PC based system that was designed as a test system for TOPAZ-II Reactor Control 
Unit (RCU) development and qualification. The simulator model has certain limitations. The 
neutronic calculations are derived from the same analysis as the TITAM; however, the 
thermohydraulic and electrical models are based upon empirical correlation using first order 
linear differential equations. The purpose behind this simplification was to allow for real- 
time computation, a requirement of RCU testing. The Real-trme Dynamic Simulator does 
not model the behavior of several subsystems, including the Electromagnetic (EM) pump, 


the batteries, and many of the gas systems. Another of its limitations 1s that the effect of 


D 


gear backlash is not included. Within these constraints, the Dynamic Simulator provides a 
reasonable approximation of TOPAZ-I behavior. Its response outside the a linear range of 


these approximations must be analyzed. 


d. Sliding Mode Control Model 


Sliding mode control is a variable structure control technique that 
uses an robust control structure. The use of high speed switching enables high accuracy, 
robust control of system response within a specially designed regime (called a sliding 
manifold). Yuri B. Shtessel developed a sliding mode control design for the TOPAZ-II 
system as part of research funded by the United States Air Force Office of Scientific 
Research (Shtessel, 1994). The mathematical modeling in the derivation of the control law 
includes a physics-based reactor kinetics model, a sophisticated model of the control drum 
drive, and a limited thermal and electrical model. The reactivity feedback model only 
includes the effect of moderator temperature changes and therefore does not account for 
the expected initially negative feedback from increasing system temperatures. The thermal 
model only tracks moderator and fuel temperatures. The electrical model assumes a 


constant cesium pressure, which does not fully account for thermionic startup conditions. 


2. Reactor Kinetics 


Al| three models use reactor kinetics formulae based upon the theory 
presented in Chapter III. (For convenience, Equations ( 6 ) and ( 7 ) are repeated below). 
The only difference is ın the treatment of Q(t), the fission equivalent source term. This ıs 
the value of initial neutron power for a reactor that has never been operated at high power. 
In the TITAM model, Q(t) is set to zero. Power initialization is accounted for elsewhere 


within TITAM programming. In both the Reactor Simulator and the Sliding Mode 


equations, Q(t) = 1.46 x IDEE y : 


iin Zu PO) pi). Vac ()* Qa) 
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where: P = reactor power (W) 


p = effective delayed neutron fraction 
p(t) = reactivity 

A, = delayed neutron fraction of ith group 

f. — decay constant of the delayed neutron precursor 

/" 2 prompt neutron lifetime 

C, = equivalent delayed neutron precursor concentration 
Q(t) = fission equivalent source power 


N = number of delayed neutron precursor groups 


di Reactivities! 


There are differences within the reactivity modeling in each simulation. All 
of the simulattons include control drum reactivity. The Sliding Mode Control Model only 
accounts for the temperature feedback reactivity of the moderator, the dominant term. The 
impact of this simplification requires analysis. A comparison of the differences among each 


model follows. 


a. Control drums 


The power regulation system sets control drum reactivity for reactor 
control. The value of the drum reactivity is initially negative since the drums contain reactor 
poison plates. Outward rotation of the drums reduces the magnitude of this negative 
reactivity, thus adding positive reactivity to the core. The control drum reactivity reported in 
TITAM is derived from personal communications with Russian researchers in. 1992 and is 
given below as Equation (10) (Paramonov and El-Genk, 1994). The source of the 
Equation (11 ), which ts used within the Reactor Simulation Model and the Sliding Mode 
Control Model, is from personal communications in 1993 (Kwok, 1994). The value of ¢, the 
initial drum reactivity, has a range of possible values due to uncertainties in the description 


of the feedback reactivity of the reactor (Shtessel, 1994). The existence of these 


! All reacuviues in this section are reported in units of 107 ($'s). 
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uncertainties makes an adaptive control scheme an appropriate choice for this type of 


system. 


D (9)94153]0:19' —357x10:* 8' +221x10° +... 


(TITAM ) (10) 


-1372x10'*8? -25x10?0 4 9 


Pp. (9) 2 689x109 - 233x107! 8* 4 328x107 8^ 4... 
(Shdng) ! (11) 


+457x10*0° - 588x1030+@ 


where: Pp,., — control drum reactivity ($) 
0 = drum angle (? ) 
¢ = initial reactivity ($) 


The resulting reactivity values for these equations are graphed in Figure 19 
below. The data reveal the maximum difference to be about ten percent of the peak value 
with similar curve shapes for each equation. The formula shown in Equation (11) will be 
used in the hybnd control model since it is the most current curve-fit in the literature. This 
is a conservative choice from a controls standpoint since the Reactor Simulator and Sliding 
Mode Control models have a lower absolute range of reactivity, which limits control 


authority. 
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Figure 19. Control Drum Reactivities. 


b. AES tor 


Moderator reactivity is an important parameter and 1s the largest 
contributor to the overall temperature reactivity. As shown in Figure 20, moderator 
reactivity is positive; therefore, an increase in moderator temperature adds positive 
feedback. The TITAM model uses a two piece fit, as given in Equation (12). Both the 
Reactor Simulator and the Sliding Mode models use the same polynomial approximation 
shown in Equation (13). The TITAM equation had a typographical error in the source 
document, report No. UNM-ISNPS-3-1994. Equation ( 12) below 1s derived from personal 
communications with the authors (Paramonov and El-Genk). Since the TITAM model has 
a slight discontinuity (see Figure 20), Equation (13) is preferred for startup control 
modeling. The Sliding Mode model has higher reactivity values, resulting in more positive 


feedback. This worst-case feedback effect justifies the choice as a cautious decision. 
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- 003027 e ERES ru e , 15:175 
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(Sliding) ( 13) 


+ 2.92 x10°(T, - 5 «4176 x10 (T, - T;) 
where: D, rar = Moderator reactivity ($) 


7,, = moderator temperature (K) 
I, = reference temperature (K) 
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Figure 20. Moderator Reactivity. 


e UO, Fuel 


The TITAM modeling uses a single function in its fuel reactivity 


formula, Equation ( 14 ), derived from 1993 research (Paramonov and El-Genk, 1994). The 


Reactor Simulator uses a two piece model, Equation ( 15). The Reactor Simulator model ts 


derived from an initial report (Gunther, 1992) and uses a linear extrapolation for values 


above 800 K. The choice of constants for this linear fit portion produced the discontinuity 


62 


shown ın Figure 21. Additionally, there 1s a large divergence between the two models in the 
calculated fuel reactivity at higher temperatures, a region of particular interest for reactor 
casualty analysis. The poor fit between the two parts of Equation ( 15) lead to the choice of 
the TITAM fuel reactivity formula for the hybrid control model. The TITAM model also 
has lower absolute values, making it a conservative choice for startup control since fuel 
temperature reactivity adds prompt negative feedback, aiding stability. Using a lower 


magnitude of fuel reactivity lowers the simulated feedback effect. 


T 
Pr 7 0432x107 (T7, - 5) 064 ee = | (14) 
(TITAM) T 


T 
Pus =12x 10T, -T)-138 |—-1|, T, «800 K (15) 
(Simulator ) h 


--10x10* -22x10 *(T, - 7), for T, 2800 K 


where: 0;,,, = fuel reactivity ($) 
T, — fuel temperature (K) 


I, = reference temperature (K) 


Note: The TOPAZ II Reactor System Real-time Dynamic Simulator 
report (Kwok, 1993) gives a value of -2.2 X 10? $/K for the slope of the fuel reactivity 
with respect to fuel temperature. The value of -2.2 X 10^ $/K is taken from the simulator 


program code directly. 
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Fuel Reactivity 
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Figure 21. Fuel Reactivity. 


d. Core Plates 


Both the TITAM and Reactor Simulator models use a linear 
approximation for the core support plate reactivity (see Equation (16) and Equation (17), 
below). The core plate reactivity 1s nearly equal to the fuel reactivity over the same range of 
temperatures. Since the time heating profile is different for each core component, the 
overall reactivity effect of core plate heating 1s reduced. This phenomenon will be analyzed 


after presentation of the other reactivity models. 


E 4 
p, =-144x10*T,-T,) T 
(TITAM ) 


= Li A = 


where: p,, = core plate reactivity ($) 
T., = core plate temperature (K) 


T, = reference temperature (K) 
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€: Reflector 

The beryllium reflector has a positive temperature coefficient of 
reactivity, as shown m Equation (18) and Equation (19). The value of this reactivity 15 
significant to control system design since reflector ejection results in the loss of this positive 


reactivity, establishing the net negative reactivity required for shutdown. 


I 18 
Pret =034 1- E (18) 


(TITAM ) ref 


Pref cosi Es (19) 
( Simulator ) T 


where: p,p = reflector reactivity ($) 
T, = reflector temperature (K) 
T, = reference temperature (K) 


£ Electrodes 


Calculation of the reactivity of the electrode is slightly more complex 
since the TITAM and Reactor Simulator models, Equations (20) and (21) respectively, 
depend upon both the emitter and collector temperatures. Figure 22 shows the calculated 


electrode reactivity profile for the first 1000 seconds of a startup. 


T T 
p, = -01 Z 1) -744 x10 E i (20) 
(TITAM ) T iP 
T. T 
p, - 85210? -426x10 | + | (21) 
Sımulator 1 I 


where: p, = electrode reactivity ($) 
T, = emitter temperature (K) 
T, = collector temperature (K) 
I, = reference temperature (K) 


Electrode Startup Reactivity 
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Figure 22. Electrode Reacuvity Dunng Startup. 


g. Reactivity Analysis 


Temperature reactivity forms the major feedback mechanism during 
startup. Consideration this reactivity is crucial to control system development. The validity 
of assumptions used within each model requires investigation. Each component heats up 
dunng startup, contributing its own reactivity value to the net core reactivity. Figure 23 
shows the net temperature reactivity feedback values calculated using the hybrid control 
model. The net reactivity line shown demonstrates the dominance of the moderator's 
temperature effect upon the system. This line is for comparison purposes only, since it give 
the net value for an isothermal core. Since the temperature differences among components 
during operation is substantial, the temperature reactivity must be examined during transient 


conditions. 
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Figure 23. Temperature Reactivity Feedback. 


The temperature changes following a step power increase from 5% 
to 95% full fission power were calculated using the TITAM simulation. The hybrid model 
reactivity equations simulates the reactivity behavior for this transient shown in Figure 24. 
This computation shows that the positive feedback effect of the moderator reactivity has a 
substantial ame delay. The net reactivity is still negative until approximately six minutes after 
the step power increase. Analysis performed for the Preliminary Nuclear Safety Assessment 
obtained similar results (Flight Safety Team, 1992). 'This provides a gauge for the fidelity of 
the hybrid reactivity equations. The Reactor Simulator model yields a comparable reactivity 
profile as long as its thermal modeling 1s adequate. Since the Sliding Mode Control Model 
only includes the moderator reactivity, its simulation would not account for the tnital 


negative feedback. The impact of this deficiency requires additional analysis. 
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Up-power Transient Reactivity Changes 
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Figure 24. Reacuvity Changes Following Up-power Transient. 


4. Thermal Model 


Among the available models, TITAM has the most complete thermal- 
hydraulic equations. The TITAM simulator results correlate well with measured data 
(Paramonov and El-Genk, 1994). The Reactor Simulator uses simple first-order exponential 
build-up equations, based upon empirical estimates of the system step-power response. The 
Reactor Simulator User’s Guide and Reference Manual does not indicate any final comparison of 
its simulation results to established data (Kwok, 1993). The Sliding Mode Control Model 
also has first order approximations for thermal modeling; however, its results were 
compared favorably with the TITAM data (Shtessel, 1994). The hybrid control model uses 
somewhat simplified verstons of the TITAM thermal equations. Also, the hybrid control 
model does not contain the additional refinements in the TITAM model that account for 
axial effects. The system temperature equations are provided roughly in order from the core 


(heat source) to the radiator (heat sink). Each symbol will be defined upon tts initial use. 
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a. Fuel Temperature 


The heat transfer relationship shown in Equation (22) is modified 
from the TITAM modeling in that it does not account for the fuel gap void. This space 1s 
between the UO, fuel pellets and the inner surface of the emitter. The fuel gap disappears as 
the fuel expands and fuses in the available space (Flight Safety Team, 1992). After the gap 1s 
filled, the heat transfer mode through the gap shifts from radiative to conductive. The 
TITAM model includes this transition stage. This effect is not included in the hybrid 


control model. 
dT, 
Mp, == PA ANT, 1) (22) 


Where M, = fuel mass (kg) 


W 
Cp, = fuel heat capacity ex) 


P = thermal power (W) 
T, = fuel temperature (K) 





W 
h, = fuel heat transfer coefficient e =) 


A, = fuel area (m?) 
T, = emitter temperature (K) 


b. Emitter Temperature 


The emitter temperature equation is derived from an energy balance 
dividing heat from the fuel between the transfer to the collector and the energy losses due 
to electron cooling and electrical lead losses. The electric power variables in Equation ( 23 ) 
are calculated with the equations described in the Thermionic Emission Model section of 


this chapter. 


dT. 
RAT M 


e m Jle (23) 
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Where M, = emitter mass (kg) 


: W 
Cp, = emitter heat capacity roa 
g . 


T, = emitter temperature (K) 
: W 
h, = emitter heat transfer coefficient pu 
m . 


A, = emitter area (m?) 

T, = emitter temperature (K) 

T, = collector temperature (K) 
P.. = electron cooling power (W) 


P,, = joule heating loss (W) 
c. Electrode Gap Temperature 
dI, 
M,Cp,— RAT TI AT, =) (24) 


Where M, = gap mass (kg) 


i W 
Cp, = gap heat capacity PS 


T, — gap temperature (K) 





W 
h, — gap heat transfer coefficient | 5 ) 
m'-K 


Á, — gap area (m?) 


d. Adjusted Gap Conductance (»,) 


The calculation of the gap heat transfer coefficient 1s somewhat 


complex it involves radiative heat transfer through cesium plasma when the system reaches 
operating temperatures. The gap contains helium gas at startup that 1s eventually replaced by 


cesium. TITAM simulates this effect with a linear exchange rate that 1s adjustable by the 
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k 
h, = +3067 T- (25) 
d, 


W 
Where A, = gap heat transfer coefficient = 
m*-K 





k ap conductivi paa 
= uctiv == 
g 7 gap ty kg- K 
o = Stephan - Boltzman's constant 
6105 ) 
| A m?.K* 


e = effective gap emissivity 
T, — gap temperature (K) 


e. Cs Vapor Conductivity (x à) 


15 _ pls 
os ST) (26) 


cds 
LH INT 
d, +114 x oe) 


cs 


Where A, = gap conductivity (wm) 
g.K 


E 


d, = gap size (m) 
T, = emitter temperature (K) 
T, = collector temperature (K) 
Fo 


S 


= cesium pressure (Torr) 


f. Cs Gap Effective Emissivity 


ef = 0142+126 x107, (27) 
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Where £^ = effective gap emissivity 
T, = gap temperature (K) 


g. Collector Temperature 


The collector temperature 1s calculated with Equation ( 28 ). This 1s 
simplified from the TITAM model by not determining the msulator and clad temperatures 
directly. The temperature differences between these elements 1s minimal due to the high 
heat transfer coefficients of the TFE materials. The collector mass and heat capacities values 
in the hybrid model are modified to account for the msulator and clad. This change 
maintains accuracy with a desirable reduction of the number of pomts in the heat transfer 


model. 


M Cp, 





T 
3 HAGO Dd P Ma 


d 
dt Jle 


Where M. = collector mass (kg) 


W 
Cp. = collector heat capacity ex) 


T, = electrode gap temperature (K) 





| W 
h, = collector heat transfer coefficient (— 3 ~) 


A, = collector area (m?) 
T, y = coolant temperature (K) 
T. = collector temperature (K) 
P.. = electron cooling power (W) 
P 


Jle 


= joule heating loss (W) 
P, = electric power (W) 


h. Coolant Bulk Temperature 





A al z I 
M yak CP Nak a -= AGC men e ue) (29) 
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Where M, — coolant mass (kg) 


W 
CP vax = coolant heat capacity = 
E . 


T, = coolant bulk temperature (K) 
W 
h, = collector heat transfer coefficient es ) 
m K 


A, = collector area (m?) 


T. = coolant temperature (K) 


m = collector flow rate (Ke) 


T., = coolant inlet temperature (K) 


5. Thermionic Emission Model 


The TITAM model includes much detailed analysis of the thermionic 
process. However, since the system startup does not include steady state control of electric 
power, certain simplifications are in order. Both the Reactor Simulator and the Sliding 
Mode Control Model use linear first order approximations. The hybnd control model 
includes simplified thermionic equations with the other power calculations required 1n its 


thermal model. 


a. Electron Cooling 


This term accounts for the heat transfer in form of electron flow 
across the interelectrode gap. As much as 20% of the thermal power 1s transferred by this 
mechanism (Flight Safety, 1992). The equation for the emitter work function, $, , 1s found 
within the Chapter III, Basic Theory. 


prega [fe a) (30) 
e e 
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- a (A 
Where / = current density (A/ 2) 
A, = emitter area (m?) 


gy, = emitter work function (eV) 


k = Boltzmann's Constant 


(1.38 x10? Jc.) 


e = fundamental unit of charge (C) 


m = collector flow rate A 


T., = coolant inlet temperature (K) 


T, = emitter temperature (K) 


b. Lead Joule Heating 


m z (JA) Ro (31) 


= f 
Where J = current density ( pa 2) 
A, = emitter area (m?) 


g. = emitter work function (eV) 
k = Boltzmann's Constant 


(1.38 x 10% Y) 


e = fundamental unit of charge (C) 


m = collector flow rate (ke) 


T, = coolant inlet temperature (K) 


T, = emitter temperature (K) 


c. Thermionic Current 


d 
z - -01X1 - 037k,P) (32) 
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Where / = thermionic current (A) 
k, = conversion efficiency 


P = reactor power (W) 


d. Thetmionic Power 
dP 
mi =-01 (2 = 037k,P) (33) 


Where P. = thermionic power (W) 
. = conversion efficiency 


P = reactor power (W) 


6. Control Drum Model 


The control drum is essential for control system design. The TITAM model 
does not include the control drum and stepper motor characteristics. The Reactor Simulator 
does not account for the hysteresis and backlash within this component. Only the Sliding 


Mode Control Model has complete equations for simulation. 


7: Hybrid Control Model 


The hybrid control model uses the equations discussed within the previous 
sections. The TITAM analysis report provided the values for most material and system 
parameters (Paramonov and El-Genk, 1994). The material properties are found within the 
Appendix of the same report. An additional approximation within the hybrid model is the 
use of linear values for many of the temperature dependent parameters such as thermal 
conductivity. The major exception to this simplification was for the thermal conductivity of 
the UO, fuel, which exhibits a large non-linear range of values over the temperatures of 
interest. For this material, the exponential approximation given in the Appendix of the 
TITAM report is incorporated with the hybrid system model. For most materials, 
temperature effects are easily accommodated by incorporating temperature dependent 


terms within the heat transfer equations. 


a. Hybrid Model Implementation 


The model is developed using Simulink®. The overall system 
diagram is shown in Figure 25. For comparison, note the correlation of this block diagram 
with the Automatic Regulating System described in Figure 17. Each system block contains 
the appropriate section of the simulation. The Matlab® function block within each group 
contain the differential equations that form the system control model. The regulating block 
calculates the appropriate control drum commands for each mode of operation. The 
control drum calculates control drum angle after receiving the stepper motor commands 
from the regulating system block. The reactivity block calculates the net reactivity based 
upon temperature feedback and control drum position. The reactor model block determines 
reactor fission power and delayed neutron precursor concentrations. The thermal model 
calculates the resulting system temperatures based upon the fission. power signal. The 


thermionics block determines the resulting electrical power parameters. 


Regulating Control Thermionic Mode 


System Drum 


Clock To Workspace 





Figure 25. Simulink& Implementation of Hybrid Control Model. 


The regulating system block presented in Figure 26 develops a 
control drum stepper motor signal. This block simulates the two basic modes of operation, 
startup and operating, and controls the stepper motor during the initial drum positioning 


during the startup sequence. 
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Figure 26. Regulating System Block. 


The regulating system block has different modes of behavior based 
for each stage within the startup. Each model includes an adjustable deadband to preclude 
system hunting. As noted earlier, limitations of the TOPAZ-II nuclear instruments result 
unreliable power level readings at low fission rates. Therefore, the initial drum positioning is 
in accordance with a programmed sequence. The drums are rotated out to 154° and then 


back in to 145°. This open loop behavior 1s accomplished with the block given in Figure 27. 
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Figure 27. Startup Drum Positioning. 
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During startup mode, the reactor power must be controlled within 
the power rate limitatons of 600 W/s, when power is less than 35 kW, and 80 W/s when 
power is above that setpoint (See Figure 16, above). This operating envelope 1s achieved 
using the Startup Power Profile block shown below, Figure 28. This block implements of 
the power rate limit by shifting to the lower power rate if power should exceed the setpoint. 
Thus, the lower limit may be set sooner than predicted 1f power overshoots the ordered 
level. Also, the higher power rate limit will be used if the power should undershoot during 
startup. Note that the overall goal of this design was to ae the envelope. There ıs 
no provision to limit power for instantaneous levels that exceed the limits. Since the Russian 
designed controller produces short periods of power rate overshoot, strict adherence is not 
essential (Shtessel, 1994). The limits themselves are based upon material concerns which are 
dependent upon the duration of heat-up rate excursions. Should further analysis indicate 
that the heat-up rates are excessive, it would be easy to incorporate such an instantaneous 


heat-up rate limiting feature within this existing design. 
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Figure 28. Startup Mode Power Profile Block. 


Figure 29 shows the construction of the Reactor Model block. The 
Reactor Simulator and Sliding Mode Models are easily investigated by replacing the 


Matlab® file for the function block with one developed for the respective model. Figure 30 
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shows simplification that may be used when the reactor is sub-prompt critical; that 1s, core 
reactivity 1s less than the value of the effective delayed neutron fraction and power 1s 
changing relatively slowly. This simplification was validated by the TITAM simulation report 
(Paramonov and El-Genk, 1994). Non fission heating in a nuclear reactor 1s primarily due to 


radioactive decay of fission products after full power operation. 
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Figure 29. Reactor Model Details. 
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Figure 30. Reactor Model - Prompt Jump Approximation Implementation. 
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b. Model Evaluation 


The following section assesses the overall effectiveness of the 
resulting model. The thermal model is validated by comparing simulated temperature 
behavior to a predetermined power history without temperature feedback. Another 
benchmark was provided by examining the data obtained from power changes performed 
upon the single TFE experiment rig. The hybrid model was compared to the TITAM 
results for a step up-power maneuver from 5% to 95% power following a system startup. A 
step power maneuver was chosen in order to isolate such non-linear effects as cesium filling 
and the ejection of the radiator. The differences between the reactivity equations for these 
models are minor; therefore, both of these simulations were examined in an open loop 
mode. This was accomplished in the TITAM simulation by simulating heating from TFE 
resistors, rather than the fission source. For the hybrid control model, the reactivity 
temperature feedback was set to zero, leaving only the control drum reactivity. As shown in 
Figure 31, there is a good match between the results of both simulations. However, the 
temperature rise following TFE loading in the hybrid model results does not agree with 
experimental data. Examination of the underlying equations reveals that this discrepancy is 
caused by a 1.6% difference in the steady state values to temperature at the 100% power 
level. Since this error is relatively small and does not significantly impact the start-up 
parameters, the simulation was not altered. As shown by the model response to the step 
power changes, the model tracks expected behavior in other respects. The fuel temperature 


response is shown; the results for the other parameters of concern are similar. 


For further analysis, the hybrid control model was replaced by the 
Reactor Simulator Model and then the Sliding Mode control model for comparison to the 
TITAM data. The Reactor Simulator Model has some discrepancies. It appears that the time 
constants chosen for this model require refinement. Figure 32 shows the equilibrium fuel 
temperatures obtained from TITAM simulations. The non-linear behavior can be expected 
since the NaK coolant flow also has electrical power dependence. The Sliding Mode 
Control Model has better correspondence with TITAM for the two temperatures it tracks. 
However, the since the Sliding Mode Control Model does not include any negative 


temperature feedback, a different analysis would be required to benchmark this simulation. 
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Thermal Model Results 
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Figure 31. Thermal Model Comparison. 
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Figure 32. Equilibnum Fuel Temperatures. 


The validity of the reactor model section and its compatibility with 
the thermal model was verified by using the net reactivity from TIT AM data as input. The 


resulting system response compared favorably, as shown in Figure 33. The power response 
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shows a noise component. These oscillations are primarily artifacts induced by the time 
intervals chosen for the integration of the differential equations. The magnitude of this 
deviation can be greatly reduced by reducing length of these time steps. However, the times 
required for calculations are greatly increased. This phenomenon is important for the 


analysis of the Sliding Mode controller, and will be considered next. 
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Figure 33. Reactor Model Evaluation. 


By restoring the temperature feedback effect, analysis of the lack of 
negative feedback within the Sliding Mode Control Model is possible. The control drum 
reactivities were used for comparison in order to track the differences between the control 
system response in both cases. Figure 34 presents the resulting control drum reactivities for 
the Hybrid model, which incorporates the negative feedback, and the Sliding Mode model, 
which does not. As can be easily seen, the magnitude of the error produced by disregarding 


negative is not significant. 
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Figure 34. Negative Feedback Comparison. 
E: CONTROL SYSTEM DESIGN 


Since the system incorporates both neutron power and current control regimes, two 
distinct control laws are required for regulating system design. There are two different 
regulation schemes available for analysis. The ongmal Russian design proposed a 
proportional feedback law for neutron power regulation and a fixed power rate change for 
payload current regulation. This programming 1s included in the Reactor Simulator code 


(Kwok). The Sliding Mode Control research includes control laws for adaptive control. 


Sliding control 1s a useful approach to systems with imprecise models or measuring. 
Such a control scheme is appropriate for the TOPAZ-II system since the feedback 
parameter, reactor power as measured by neutron current, can have a large degree of 
uncertainty (on the order of 50%). Since these detectors cannot be calibrated after launch, 
their accuracy over time 1s uncertain. Sliding control 1s a robust control scheme that applies 
an optimization process upon the error function. This error function 1s modeled as a series 
of nth order differential equations. The optimization process is obtained by minimizing the 


integral of time multiplied by the absolute error. For neutron and current regulation, the 
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sliding surfaces given in Equation ( 34 ) were derived from the hybrid equations developed 


above. 


c, 2 5" 012599 «0327 20 


(34 
O, =€"+2.68£"+1944g'+0.639£=0 i 


where: 7=P-P area 
= power tracking error 


eT. 


— current tracking error 


Both switching control laws are synthesized from U = V,,,,*sign(6), where U is the 
motor control signal which 1s a step signal with values between +27 to -27 V. Since these 
control laws involve higher order derivatives and the sign operator, a non-linear function, 
the system control must be carefully designed so as to prevent chattering. This behavior is 
kept in check by smoothing the control law, which improves system robustness. Within 
Matlab®, this is readily implemented with a saturate block. Other options include increasing 
the system sampling frequency, which increases the resulting control response calculations. 
However, this additional computational overhead 1s not desirable in view of the requirement 
to test the control system with a simple microprocessor. Figure 35 gives the system used 
within the Hybrid model, where hunting has been reduced with the use of a saturation 
block. The sliding mode control report did not include the characteristics of the method 
used in his analysis (Shtessel, 1994). The trade-off in adjusting the control law is that some 
tracking precision and bandwidth 1s lost. (Slotine and Li, 1991) 
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Figure 35. Sliding Mode Control Implementation. 
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D. DESIGN EVALUATION 


For both control schemes, a standard startup was evaluated with the closed system 
temperature feedback enabled. The results for the proportional control and sliding mode 
control schemes are shown in Figure 36. Examination of both curves show that the sliding 
mode control scheme has a faster response to such an input. The proportional control 
exhibits more noticeable lag, but does not show significant overshoot or undershoot. The 
instantaneous power rate exceeds the startup limit for a few short periods. Both control 
laws track the ordered reactor power with good precision. The sliding control proved to be 
very sensitive to the range selected for the saturation. Without proper values, the random 
noise like signal produced by the quantization of the simulator was enough to produce 
chattering. Further evaluation revealed that a random noise signal of ~5% of power signal 


was enough to induce severe oscillations. 
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Figure 36. Control System Evaluation. 
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E SCENARIO ANALYSIS 


In this section, performance of the hybrid control system design 1s considered in the 
case of different casualties and startup conditions. This step 1s critical in design acceptance 


in order to meet the goal of reducing risk to as low as reasonably achievable. 


1. Reactivity Addition Casualty 


Addition of large amounts of positive reactivity will cause a power 
excursions. The worst case casualty considered was a failure of the control system, leading 
to control drum withdrawal to the limit. Figure 37 shows the result of such a power 
excursion upon the hybrid simulation model. This casualty was the first analyzed since it was 
also considered in the Preliminary Nuclear Safety Assessment, which provided for easy 


comparison between the two models. 
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Figure 37. Control System Failure. 
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The dramatic effects revealed by the data shown in Figure 37 reveal that fuel 
temperatures most quickly rise to meet the limits. Therefore, the behavior of fuel 
temperature provides a good indication of the adequacy of control system safety response. 
With temperature limiting safety feature in effect, thermal limits are not reached for this 
type of casualty, as shown in Figure 38. Since the complete simultaneous failure of both the 
control system and the associated safety temperature limiting is unlikely, no further 
protection of the system 1s required. Sufficient protection for such a casualty would be easily 
accommodated through the use of redundant safety circuits, a standard practice within 


reactor control system design. 
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Figure 38. Temperature Limiting Response. 
2 Loss of Coolant Accident 


In thıs casualty, the reactor response is considered after a coolant system leak. For ease of 
analysis, an exponential loss of coolant flow was assumed followed by a complete loss of 
coolant flow. An actual loss of coolant would be more severe, since the remaining NaK 


would provide some thermal heat capacity to absorb heat. For this case, in accordance with 
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the actual loss of coolant casualty experience during the Russian experiments, it 15 assumed 
that the coolant limiter would not be effective in initiating safety system response. As 
shown in Figure 39, limiting core temperatures are reached, making this casualty of 


particular concern for overall system evaluation. 
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Figure 39. Loss of Coolant Casualty Analysis. 

Since the temperature limiter is the only operational safety feature originally 
planned for the TOPAZ-II system, system modification would be required. For a loss of 
coolant casualty, there are several options to initiate reflector ejection. The NaK coolant 
pressure bellows could provide an input upon loss of pressure. However, mechanical 
bellows are not as reliable as purely electronic means. Another option would be to track 
differential coolant temperature. An anomaly with coolant flow would be detectable 
through this means with some time delay. The fastest acting option would be to use means 
that could directly detect the presence of the coolant. Since the coolant is a liquid metal, 
resistance detection seems a reasonable solution. As shown in Figure 40 below, fuel 


temperatures can be kept below limits 1f the reflector eJection is initiated in. time. 
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Coolant Loss with Reflector Ejection 
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Figure 40. Coolant Loss with Protective Action. 


3. Cesium Pressure Oscillations 


The presence of cold spots within the Cs system will cause undesired 
oscillanons in the output current. The main consequence with for the control system is that 
the system should not attempt to maintain. fission. power during such a swing in current 
output. Such a response would cause additional thermal power changes, which would be 
likely to cause further oscillations. The option examined within the hybrid simulation was to 
also track differential coolant temperatures in order to provide a backup means of 
maintaining the constant power output. As shown in Figure 41, such a strategy does prevent 
the control system from following the output current instability. Additional options would 
include delaying current following mode until current oscillations had ceased. Experience 
with the TOPAZ test facility showed that the system would stabilize without operator 
action. However, since the presence of such current swings was a newly observed 
phenomenon, it can not be reliably predicted. Therefore, delaying onset of current tracking 


mode would require careful control by the ground station. 
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Figure 41. Cesium Pressure Induced Oscillations. 
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VIII. CONCLUSIONS 


Analysis of the startup the TOPAZ-II nuclear space reactor indicates that design of 
a control system to meet required performance standards 1s realizable. In particular, sliding 
mode control is well suited for use within such a system. Regardless of the control method 
chosen, careful attention must be give to the simplificanons and approximations used. The 
startup system behavior includes many non-linear elements. The effectiveness of a hybrid 
simulation that can incorporate these elements 1s shown. The assumptions contained within 
the dynamic simulator model require further refinement to become a more reliable test 


mechanism for evaluating startup control programming. 


The sliding mode control method of power regulation 1s shown to be an efficient 
means of establishing the desired reactor power profiles. The sensitivity of the controller to 
non-linear measurement is a key concern. The noise sensitivity 1s an important 
consideration especially in view of the high radiation environment present both in space and 


around an operating nuclear core. 


The physical robustness of the system was demonstrated following the loss of 
coolant casualty. The preliminary analysis presented here provides groundwork for 
evaluating control system changes with respect to reducing risks to as low as reasonably 


achievable. 


ue 





IX. RECOMMENDATIONS 


The Reactor Simulator model could be improved by incorporaung some of the 
elements effectively demonstrated within the Hybrid Model presented here. The Matlab® 
implementation of the hybnd simulation provided real-trme results when compiled. This 


indicates that the requirement for real-time evaluation can be maintained. 


The noise characteristics of the system require further evaluation to ensure the 
proper selection of controller method. This is particularly true in the low power regimes, 
where the detectors will be providing a noisier signal. The installation of low range neutron 


detectors would provide a means of reducing this concern substantially. 


The initiation of reactor shut down following a loss of coolant casualty requires 
system modification. The use of a resistance based coolant detection system is proposed. 
Additionally, the monitoring of coolant differential temperatures would provide a secondary 


means of monitoring system thermal power output. 





APPENDIX. SIMULINK& PROGRAM CODE 


The Sımulınk® code files used in the closed loop hybnd simulation 1s presented 


here for reference. 



































function [ret,x0,str,ts,xts]=closed(t,x,u,flag); 
SCLOSED is the M-file description of the SIMULINK system 
named CLOSED. 
% The block-diagram can be displayed by typing: CLOSED. 
% 
% SYS=CLOSED(T,X,U,FLAG) returns depending on FLAG certain 
% system values given time point, T, current state vector, 
X, 
% and input vector, U. 
% FLAG is used to indicate the type of output to be returned 
in SYS. 
% 
% Setting FLAG=1 causes CLOSED to return state derivatives, 
FLAG=2 
% discrete states, FLAG=3 system outputs and FLAG=4 next 
sample 
% time. For more information and other options see SFUNC. 
% 
% Calling CLOSED with a FLAG of zero: 
% [SIZES]=CLOSED([],[],[],0), returns a vector, SIZES, which 
% contains the sizes of the state vector and other 
parameters. 
SIZES(1) number of states 
2) number of discrete states 
SIZES(3) number of outputs 
4) number of inputs 
5) number of roots (currently unsupported) 
6) direct feedthrough flag 
7) number of sample times 


de ge ge ge ge ge ge ge ge 
o 
HH 
N 
m 
ao 


For the definition of other parameters in SIZES, see 
SFUNC. 
% See also, TRIM, LINMOD, LINSIM, EULER, RK23, RK45, ADAMS, 
GEAR. 


% Note: This M-file is only used for saving graphical 
information; 
% after the model is loaded into memory an internal model 


oF 


% representation is used. 


$; the system will take on the name of this mfile: 
Sys - mfilename; 
new system(sys) 
simver(1.3) 
if (0 == (nargin + nargout)) 
set param(sys,'Location',[99,158,739,491]) 
open system(sys) 
end; 
set_param(sys,'algorithm', 'RK-45') 
set_param(sys,'Start time', '0.0') 
set_param(sys,'Stop time', '3000') 
set_param(sys,'Min step size', '0.001') 
set_param(sys,'Max step size', '.2') 
set_param(sys,'Relative error','1e-3') 


set_param(sys,'Return vars', '') 
set param(sys,'AssignWideVectorLines','on'); 


[add block('built-in/Note',[sys, '/' ,'Reactivity']) 
set param([sys,'/','Reactivity'],... 
‘position 265,70,270 751) 


add block('built-in/Note',[Sys, /', 'Power']) 
set param([sys,'/','Power'],... 
'position',[395,70,400,75]) 


% Subsystem 'Core Reactivity'. 


new_system([sys,'/','Core Reactivity']) 
set param([sys,'/','Core 
Reactivity'],'Location',[0,59,274,252]) 


add block('built-in/Inport',[sys,'/', 'Core Reactivity/x']) 
set param([Ssys,'/','Core Reactivity/x'],... 
position 655,55 ,85,75]) 


add block('built-in/S-Function',[sys,'/',['Core 
Reactivity/S-function',13,'M-file 
plots 3 Ines 195." T) 
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set param([sys,'/',['Core Reactivity/S-function' ,13, 'M-file 
which plots',13,'lines MID 

"function name','sfunyst',... 

'parameters','ax, color, npts, dt',... 

'position',[130,55,180,75]) 
add_line([sys,'/','Core Reactivity'],[90,65;125,65]) 
set_param([sys,'/','Core Reactivity'],... 

'Mask 
Display','plot(0,0,100,100,[83,76,63,52,42,38,28,16,11,84,11 
,11,11,90,90,11],[75,58,47,54,72,80,84,74,65,65,65,90,40,40, 
85790]) ',... 

'Mask Type','Storage scope.') 
set param([sys,'/','Core Reactivity'],... 

‘Mask  Dialogue','Storage scope using MATLAB graph 
window.\nEnter plotting ranges and line type.|Initial Time 
Range:|Initial y-min:|Initial y-max:|Storage pts.:|Line type 
(rgbw-.:xo):') 
set_param([sys,'/','Core Reactivity'],... 

'Mask Translate','npts = @4; color = @5; ax = [0, @1, @2, 
03]; dt=-1;') 
set param([sys,'/','Core Reactivity'],... 

'Mask Help','This block uses a MATLAB figure window to 
plot the input signal. The graph limits are automatically 
scaled to the min and max values of the signal stored in the 
scope''s signal buffer. Line type must be in quotes. See the 
M-file sfunyst.m.') 
set param([sys,'/','Core Reactivity'],... 

' Mask 
Entries','5\/-10\/10\/200\/''y-/g--/c-./w:/m*/ro/b+''\/") 


$; Finished composite block 'Core Reactivity'. 


set param([sys,'/','Core Reactivity'],... 
'position',[280,165,310,205]) 


% Subsystem 'Reactor Model '. 





new_system([sys,'/','Reactor Model']) 
set param([sys,'/','Reactor 
fedele], 'Location',[10,79,751,471]) 
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add block('built-in/Note',[sys,'/', ‘Reactor Model/Reactor 
Kinetics Model']) 
set _param([sys,'/','Reactor Model/Reactor Kinetics 
Model'],... 

'Font Weight','bold',... 

'Font Underline',1,... 

Font o1Ze 26, .. 

"BoS1t10n, 5,1320;5,325,10]) 


add block('built-in/Inport',[sys,'/', Reactor 

Model/Reactivity']) 

set param([sys,'/','Reactor Model/Reactivity'],... 
'position',[55,165,75,185]) 


add_block('built-in/Mux',[sys,'/',['Reactor 
Model/State',13,'Variables']]) 
set param([sys, '/',['Reactor 
Model/State',13,'Variables']],... 
PIBDDUtS ..2',... 
'position',[160,136,190,169]) 


add block('built-in/Integrator',[sys, '/', Reactor 
Model/Integrator']) 
set param([sys,'/','Reactor Model/Integrator'],... 
EInyrtial "(Come en 
‘position’ ,[450,140,480,170]) 


add block('built-in/MATLAB Fen',[sys,'/',[ ‘Reactor 
Model/Reactor',13,'Calculations']]) 
set param([sys, '/',['Reactor 
Model/Reactor',13,'Calculations']],... 

‘Drop Shadow' ,4,... 

SMATEABZEEN | CX pa s 

OU PUt Width! EN 

'position',[265,136,340,174]) 


add_block('built-in/Note',[sys,'/','Reactor Model/Delayed 
Neutron Precursor Concentrations']) 
set_param([sys,'/','Reactor Model/Delayed Neutron Precursor 
Concentrations'],... 

“pesition. ,[285,95,290, 100] 


add block('built-in/Note',[sys,'/', 'Reactor Model/Prompt 
Jump Approximation Calculations']) 
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set param([sys,'/','Reactor Model/Prompt Jump Approximation 
Calculations db s. 
LEON size dum... 









dposition',[315,50,320-55]0 





add block('built-in/MATLAB 
Model/Neutron',13,'Power']]) 

set_param([sys,'/',['Reactor Model/Neutron',13,'Power']],... 
'Drop Shadow',4,... 

'MATLAB Fcn','P pja',... 
Euu:put Width','1',... 
'position',[480,201,555,239]) 


Fen',[sys,'/',['Reactor 











add block('built-in/Outport',[sys, '/', ['Reactor 
Model/Reactor',13,'Power']]) 

set param([Sys,'/',['Reactor Model/Reactor',13,'Power']],... 
sposition’ ,(685,205, 705, 225] ) 







add block('built-in/Sum',[sys,'/','Reactor Model/Sum']) 
set param([Sys,'/','Reactor Model/Sum'],... 
'position',[620,205,640,225]) 









add block('built-in/Constant',[sys, /', ['Reactor 

Model/Non-fission',13,'Heating']]) 

set param([sys,'/',['Reactor 

Model/Non-fission',13,'Heating']],... 
parye','Qo',... 

'position',[525,127,565,153]) 
add_line([sys,'/','Reactor 
Model'],[80,175;135,175;135,160;155,160]) 
add line([SyS,'/','Reactor 
Model'],[485,155;495,155;495,130;140,130;140,145;155,145]) 
add line([Ssys,'/','Reactor Model'],[195,155;260,155]) 
add line([sys,'/','Reactor Model'],[345,155;445,155]) 
add line([sys,'/','Reactor Model'],[645,215;680,215]) 
add line([sys,'/','Reactor Model'],[560,220;615,220]) 
add line([syS,'/','Reactor 
Model'],[570,140;600,140;600,210;615,210]) 
add line([sys,'/','Reactor 
meget | ,(225,155;225,220;475,220] ) 

















101 


*; Finished composite block 'Reactor Model'. 


set param([sys, ' /','Reactor Model'],... 
'position',[325,69,355,121]) 


% Subsystem ['Control',13,'Drum',13,'Position']. 


new_system([sys,'/',['Control',13,'Drum',13,'Position']]) 
set param([sys,'/',['Control',13,'Drum',13,'Position']],'Loc 
ation',[27,43,780,497]) ; 


add_block('built-in/Note',[sys,'/',['Control' ,13,'Drum' ,13, ' 
Position/Drum Model']]) 
set param([sys,'/',['Contro1',13, Drum',13, 'Position/Drum 
Model']],... 

Sposa ron [1505525155510] 


add_block('built-in/Inport',[sys,'/',['Control',13,'Drum',13 
, Position/Reactor',13,'Power']]) 
set param([sys,'/',['Control',13,'Drum',13,'Position/Reactor 
wis, Power ll... 

'position',[30,40,50,60]) 


add _block('built-in/Relay',[sys,'/',['Control',13,'Drum',13, 
'Position/5 kW',13,'Detection']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/5 
kW',13,'Detection']],... 

'On switch value','5000',... 

'Off switch value','-5000',... 

'position',[115,188,145,212]) 


% Subsystem ["Control',13,'Drum' ,13,'Position/In ti 
Out’ ,13,'Switcn’ |. 


new_system([sys,'/',['Control',13,'Drum',13,'Position/In // 
] 


Out 3 13, switch’ | }) 
seteoparam([sys, /' ,.| Contool l3. Doum 313, Position, in fy 
Out' ,13, 'Switch']],'Location' ,[159,417, 467,586] ) 


ada Dblock{ ‘built-in/Outport’, (sys, / ,{ Control” 13, Drumaes 
37,595 LE Tonne / Out as swetch/ out 1"]]) 
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set param([sys,'/',['Control' ‚13, Positron; in l 
Out',13,'Switech/out 1i m 


'position',[265,70,285,90]) 


,13,'Drum' 








add_block('built-in/Inport',[sys,'/',['Control',13,'Drum',13 
,' Position/In // Out',13,'Switch/in 1']]) 

set param([sys,'/',['Control',13,'Drum',13, 'Position/In " 
Buts 13, Switch/in 1 |] mm 

iposition',[35,30,55,50]D 













add block('built-in/Relational 
Operator',[sys,'/',['Control',13,'Drum',13,'Position/In Gh 
Out',13,'Switch/Relational',13,'Operator']]) 
Mampanen(Isyss Js [Controle D oUm: bool tren an 127 
Out',13,'Switch/Relational',13,'Operator']],... 

‘operator’, ENE 
'position' 1140, ‚32, 170,63]) 










ESdSblock('built-in/Sum',[sys, ' /',['Control',138, 'Drum ,13, P 
osition/In // Out',13,'Switch/Sum']]) 

BENED cons ME // 
BUTS, 13, M ee 
‘inputs’, ota 

‘position: ‚[215, 64, 235,91]) 














add block('built-in/Relational 
BEDonstor',[sys,'/',['Control',13,'Drum', 13, 'Position/In if 
Out' ,13, 'Switch/Relational' ,13, 'Operator1']]) 
MEUM > 2s | Control", 13, Drum", 13, Position/In Ph 
Out' ,13, 'Switch/Relational' ,13,'Operatori1']],... 

‘operator’, OP 
'position' 1140, 92, 170,123] ) 





add_block('built-in/Constant',[sys,'/',['Control',13,'Drum', 
13,'Position/In // Out',13,'Switch/Constant']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/In // 
Ut 13, Switch/Constant']],... 

US 

'position',[65,105,85,125]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/In ti 
Out' ,13, 'Switch']],[60,40;135,40]) 
add _line([sys,'/',['Control',13,'Drum',13,'Position/In / / 
ue, 13, switch']],[95,40;95,100; 135,100] ) 
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add line([Sys,'/',['Control',13,'Drum', 13, 'Position/In 
Out',13,'Switch']],[90,115;135,115]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/In // 
Das waitch 17, 11107115:110,55;135,55]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/In Ij 
poma switch], [1715 1/05 185,110:185,85;210,85]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/In ia 
Oui 3 5 Swxbchogs175,50:185,50:185,70;210, 7/01) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/In // 
Out' ,13, 'Switch']],[240,80;260,80] ) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/In iy 
Quis > SECHE 

' Mask 
Display','plot(-50,-50,50,50,[-50,50],[0,0],[0,0],[-50,50],[ 
-40,0],[-30,-30],[0,40],[30,30])',... 

'Mask Type','Sign',... 

‘Mask Dialogue','y = sign(x)') 
set param([sys,'/',['Control',13,'Drum', 13, 'Position/In 2 
QUE Swi teh |. 

'Mask Help','Sign Function:\n\t\t\ty = 1 if x > O\n\t\t\ty 
= 0 if x = O\n\t\t\ty = -1 if x < 0') 


% Finished composite block 
[ Control',13, Drum’ ,13, 'Position/In // Out',; 13, Switch jm 













set param([sys,'/',['Control',13,'Drum',13, 'Position/In ul 
Out 713; Swatch] 
‘position’ , [475,187 ,505,213] ) 


% Subsystem 
['Control' ,13, 'Drum' ,13, 'Position/Drum' ,13, 'Response' ,13, 'Mo 
del']. 





new_system([sys,'/',['Control',13,'Drum',13,'Position/Drum', 
13,'Response',13,'Model']]) 

set param([sys,'/',['Control',13,'Drum', 13, 'Position/Drum',1 
3,'Response',13,'Model']],'Location',[94,85,665,353]) 
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add block('built-in/Productibsys, 77,17 Control 137 D run 
3,'Position/Drum',13,'Response',13,'Model/Product']]) 

set param([sys,'/',['Contro1',13, 'Drum' , 13, 'Position/Drum',1 
3, 'Response',13,'Model/Product']],... 
'position',[230,143,260,167]) 










add block('built-in/Integrator' ,[sys,'/’,[’Control (13, Drum 
',13,'Position/Drum' ,13, 'Response' ,13, 'Model/Integrator' ]]) 
set_param([sys,'/',['Control' ,13, 'Drum' ,13,'Position/Drum' , 1 
3,'Response',13,'Model/Integrator']],... 
'position',[315,142,345,168]) 



























Bagdesbiock('built-in/Outport',[sys, / ' ,[wContro?',139, Drum 1 
3,'Position/Drum' ,13, 'Response' ,13, 'Model/Drum' ,13, 'Position 
B 

set param([sys,'/',['Contro1',13, 'Drum', 13, 'Position/Drum',1 
3,'Response',13,'Model/Drum',13, 'Position']],... 
'position',[440,145,460,165]) 


add block('built-in/Constant',[sys,'/',['Contro1',13, 'Drum', 
13,'Position/Drum',13,'Response',13,'Model/Drum Speed']]) 
EEBEDaram([sys, /',['Control',13,'Drum' , 18, "Position/Drum',1 
3, Response',13,'Model/Drum Speed']],... 

dordentation',2,... 
E21ue','1.4',... 
iposrtion',[240,77,275,103]) 


EBddEbIOSK(' built-in/Inport',[sys,'/',['Control',13, 'Drum' , 13 
, Position/Drum',13,'Response' ,13, 'Model/ In [1 
BUt 13, 'Switch']]) 
Eecgpapam([sys, ' /' ,['Control',19, ' Drum' ,13, 'Position/Drum',1 
3,'Response',13,'Model/In // Out',13, Switch ]],... 

sees tion’ ,[(80,150, 100% 170i) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Drum',13 
,s Response',13,'Model']],[235,90;120,90;120,150;225,150]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Drum',13 
mEsponse',13,'Model']],[265,155;310,155]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Drum',13 
, Response',13,'Mode1']],[350,155;435,155]) 
Egeme([sys,'/',['Control',13,'Drum', 13, 'Posrtion/Drum*, 13 
, Response',13,'Model']],[105,160;225,160]) 
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% Finished composite block 
['Control' ,13,'Drum' ,13, 'Position/Drum' ,13, 'Response' ,13,'Mo 
del']. 


set_param([sys,'/',['Control' ,13, ‘Drum' ,13,'Position/Drum' ,1 
3,'Response',13,'Model']],... 
"posrtiom4535 175.565,225]) 


add _block('built-in/Saturation',[sys,'/',['Control',13,'Drum 
',13,'Position/Drum',13,'Limits']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Drum',1 
Sintesi] 

eN c. 

‘Upper Limit','18000',... 

‘position’ ,[610,188,640, 212] ) 


add block{  bui@t-in/Outport’ ,[sys,’/",[ “Control”, 13, Drums: 
3,'Position/Drum Angle']]) 
set param([sys, '/',['Control',13,'Drum', 13, 'Position/Drum 
Angle']],... 

'position',[700,190,720,210]) 


% Subsystem 
['Control',13,'Drum',13,'Position/Drum',13,'Angle']. 


new system([sys,'/',['Control',13,'Drum',13, 'Position/Drum', 
13, 'Angle']]) 

set param([sys,'/',['Control',13,'Drum',13, 'Position/Drum',1 
3,'Angle']],'Location',[0,59,274,252]) 


add block( built-in/Inport' ,[sys, /',['"Control',13, Drum ie 
, Position/Drum',13,'Angle/x']]) 
set_param([sys, '/',['Control", 13, 'Drum' , 13) ' Position/Drum 5M 
3, 'Angle/x']],... 

'position',[65,55,85,75]) 


add block('built-in/S-Function',[sys, /',['Control',13,'Drum 
',13,'Position/Drum',13,'Angle/S-function',13,'M-file which 
plots 4.13, lines 218.271 
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set param([sys,'/',['Control',13,'Drum', 13, 'Position/Drum',1 
3,'Angle/S-function',13,'M-file which 
plots' ,13, Lines ..19, Jj PE 

'function name','sfunyst',... 

'parameters','ax, color, npts, dt',... 

iposition', [130,55,180 za) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Drum',13 
,'Angle']],[90,65;125,65]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Drum',1 
E Angle']],... 

Mask 
pisplay','plot(0,0,100,100,[83,76,63,52,42,38,28,16,11,84,11 
,11,11,90,90,11],[75,58,47,54,72,80,84,74,65,65,65,90, 40,40, 
90,90])',... 

'Mask Type','Storage scope. ') 





Set param([sys,'/',['Control',13,'Drum',13,'Position/Drum',1 
3,'Angle']],... 

‘Mask Dialogue','Storage scope using MATLAB graph 
window.\nEnter plotting ranges and line type.|Initial Time 
Range: |Initial y-min:|Initial y-max:|Storage pts.:|Line type 
(rgbw-.:xo):') 





EEcMp3pam([syS,'/',['Control',13, ' Drum , 13, 'Position/Drum';1 
er angie')]],... 

‘Mask Translate','npts = 64; color = @5; ax = [0, @1, @2, 
Q3]; dt=-1;') 

Eparam(lsys,' /',['Control',13;" Drum" 943, “Position/Drum ,1 
Exngle']],... 

'Mask Help','This block uses a MATLAB figure window to 
plot the input signal. The graph limits are automatically 
scaled to the min and max values of the signal stored in the 
scope''s signal buffer. Line type must be in quotes. See the 
M-file sfunyst.m.') 

ESwapegpam([sys, '/',['Control',13,'Drum', 13, 'Position/Drum',1 
Sm Angie ]],... 

'Mask 

emeptes  ,'5\/-10\/10\/200\/''y-/g--/c-./w:/m*/ro/p+' '\/') 
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% Finished composite block 
Control. .13, Drum’ ,13;,'Position/Drum',13, 'Angle']). 


set_param([sys,'/',['Control',13,'Drum',13,'Position/Drum',1 
3 Angle“ J]; 
'position',[690,65,720,105]) 


% Subsystem 
[Control 13, Drum 13, Position/Startup',13,'Drum', 13, POS 
it LOM Ge lhe . 


new_system([sys,'/',['Control' ,13, 'Drum' ,13, 'Position/Startu 
pgs Drum, 13, Positioning ]]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Startup 
',13,'Drum',13,'Positioning']],'Location',[24,152,771,460]) 


add_block('built-in/Qutport',[sys,'/'",["Control',13, Drúmii 
3,'Position/Startup',13,'Drum',13,'Positioning/Drum',13,'Mot 
ones. switch |") 

set_param( [sys,'/'„[.'Control',13,'Drum' , 13,2 Position/Starzeg 
-,13,'Drum’,18, 'Positioning/Drum' ‚13, "Motor" ‚13, 'Switeh Ye 


POSITION [6555 1407675, 160]) 


add block('built-in/Dead 
Zone',[sys, /',['Control',13,'Drum', 13, 'Position/Startup PEE: 
, Drum', 13, 'Positioning/Drum',13,'Position',13,' Dead 
Band']]) 
set param([sys, ' /',['Control',13,'Drum', 13, 'Position/Star tp 
',13,'Drum' ,13, 'Positioning/Drum' ,13,'Position' ,13, 'Dead 
Band']],... 

'Lower_value','-0.01',... 

'Uppere value’ ;s0.01',..= 

'position',[560,137,590,163]) 


add_block( 'built-in/Sum',[sys,'/',['Control' ‚13, 'Drum‘ , 1388 
osition/Startup',13,'Drum',13,'Positioning/Position',13, ‘Em 


or']]) 
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set param([sys, P | Gonup Drum 1, o it on S 
13. DUNAS Position LE 
“inputs' x 

‘position’ ‚[485, 140, 505,160] ) 





add block('built-in/Switch',[sys, /',['Contro1',13, Drum' ,13 
,' Position/Startup',13,'Drum', 13, 'Positioning/Drum',13,'Post 
pon' ,13, 'Switch' ]]) 

set param([sys,'/',['Control',13,'Drum',13, 'Position/Startup 
,13,'Drum',13,'Positioning/Drum',13,'Postion',13,'Switch']] 
E. 

'Threshold','1', 

'position',[340,94,370,126]) 






















add_block('built-in/Inport',[sys,'/',['Control',13,'Drum',13 
, Position/Startup',13,'Drum',13,'Positioning/Drum',13,'Angl 
e']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Startup 
Eis Drum',13,'Positioning/Drum',13,'AÀngle']],... 
gposition',[25,25,45,45]) 























add _block('built-in/Constant',[sys,'/',['Control',13,'Drum', 
13, 'Position/Startup',13,'Drum',13,'Positioning/Position',13 
‚'#1']]) 
set param([sys,'/',['Control',13,'Drum',13, 'Position/Startup 
E bDrum',13,'Positioning/Position ,13,'81']|]$—-.. 
'Value','154',... 
ipesition',[225,165,275,185]) 





addgblock('built-in/Constant',[sys, '/',['Control',13, Drum, 
IS Position/Startup',13, 'Drum', 13, 'Positioning/Position',13 
Exec) ) 
3 
set param([sys,'/',['Control',13,'Drum', 13, 'Position/Startup 
Er Drum',13, Positioning/Position',13,'482']],... 
'Value','145', 
position’ ,(225,35,275,55]) 





add_block('built-in/Relay',[sys,'/',['Control',13,'Drum',13, 
'Position/Startup',13,'Drum',13,'Positioning/154 
, 13, 'Degree',13, 'Relay']]) 
gEeEEparam([sys,' /',['Controk',13,'Drum",; 13, 'Position/Startup 
Bis Drum’ ,13, 'Positioning/ 154 
wl, Degree’ ,13, Relay']]),... 
'On switch value','154', 
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'Qff switch value','O',... 

sOnmoutput Value’, °10'5... 

position .[155,98,185,122 ] ) 
agdmline(isys, /',['Control’,13, 'Drum',13,'Position/Startup' 
, 13, 'Drum',13, 'Positioning']],[280,45;300,45;300,100;335,100 
]) 


add line qisys,-/ >| Control 413, ‘Drum' ,13, 'Position/Startup! 
‚13, Drum 5e. Posrtrionsmg' ]],[280,175;300,175; 300,120; 335,1 
20]) 

add line([Sys,'/',['Control',13,'Drum',13, 'Position/Startup' 
,13,'Drum' ,13,'Positioning']],[190,110;335,110]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Startup' 
,13,'Drum',13,'Positioning']],[375,110;420,110;420,145;480,1 
45]) 


add_line([sys,'/',['Control',13,'Drum',13,'Position/Startup' 
"13 Drum , 13, Positiemina')]].,[50,35565,35;65, 110150 AO» 


add_line([sys,'/',['Control',13,'Drum',13,'Position/Startup' 
,13,'Drum',13,'Positioning')]),[595,150;650,150]) 

add _line([sys,'/',['Control',13,'Drum',13,'Position/Startup! 
19, ' Drum',13, Positsening']],(510,150;555, 150]) 

add _line([sys,'/',['Control',13,'Drum',13,'Position/Startops 
,13,'Drum',13,'Positioning']],[65,110;65,260;420,260;420,155 
:480,155] ) 


% Finished composite block 
[ Control’ ,13,'Drum' ,13, 'Position/Startup ,13, 'Drum’,13, Fae 
itiománg"]. 


set_param([sys,'/',['Control',13,'Drum',13,'Position/Startup 
ODUM Ts Positioning l| EE 
“position (1185250 215 3001) 


add_block('built-in/Switch',[sys,'/',['Control' ,13,'Drum' ,13 
, Position/Startup',13,'Mode', 13, 'Switch'1]) 
set param([sys,'/',['Control',13,'Drum', 13, 'Position/Startup 
Eo Mode ^13 "Switeh!]] = 

anumeshold'o 1... 
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'position',[385,184,415,216]) 


% Subsystem 
['Control' ,13, 'Drum' ,13, 'Position/Power' ,13, 'Rate' ,13, “Posit 
Bonang ]. 


new system([sys,'/',['Control',13,'Drum', 13, 'Position/Power' 
,13,'Rate',13,'Positioning']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13,'Rate',13,'Positioning']],'Location',[40,48,705,323]) 


% Subsystem 
['Control',13,'Drum',13,'Position/Power',13,'Rate',13, 'Posit 
ioning/Startup',13,'Power',13, 'Profile']. 


new_system([sys,'/',['Control',13,'Drum',13,'Position/Power' 
,13, 'Rate' ,13, 'Positioning/Startup' ,13, 'Power' ,13, 'Profile' ] 
]) 
BEeBEparam([sys,'/',['Control',13,'Drum', 13, 'Position/Power', 
13,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile']] 
, Location',[27,43,780,476]) 


aqd bDlock('built-in/Inport',[sys,'/',['Control',13,'Drum',13 
, Position/Power',13,'Rate',13,'Positioning/Startup',13,'Pow 
er',13,'Profile/Reactor',13,'Power']]) 

set param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13, 'Rate',13, 'Positioning/Startup',13,'Power',13,'Profile/Re 
Ere 13, Power']],... 


position", [30,10,50,30]) 


agdEbIock('built-in/Product',[sys, '/',['Contro1*,13, 'Drum' , 1 
3, Position/Power',13,'Rate',13, 'Positioning/Startup',13, 'Po 
wer',13,'Profile/Product']]) 
sew param([sys,'/',['Control' ,13, ‘Drum’ ,13, 'Position/Power' , 
harte” ‚13, 'Positioning/Startup' ‚13, 'Power',13,'Profile/Pr 
Dpoguct ]]1,... 

'position',[390,209,420,251]) 
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add block('built-in/Inport',[sys, /',['Control',13,'Drum' ,13 
, Position/Power',13,'Rate',13, 'Positioning/Startup',13, 'Pow 
er',13,'Profile/Power',13,'Rate',13,'Limit',13,'Enabled']]) 
set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13, 'Rate' ,13, 'Positioning/Startup' ,13, 'Power' ,13, 'Profile/Po 
vena Have .13, Limit ’,16, Enabled']],... 

A A 

position” ,[330,230,350,250]) 





























add_block('built-in/Integrator',[sys,'/',['Control',13,'Drum 
',13,'Position/Power',13,'Rate',13,'Positioning/Startup',13, 
'Power',13,'Profile/Ordered',13,'Power',13,'Calculation']]) 
set param([sys,'/',['Control',13,'Drum', 13, 'Position/Power', 
13,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Or 
dered',13,'Power',13,'Calculation']],... 

initial. 5000" 1... 
'position',[465,220,505,240]) 





add _block('built-in/Constant',[sys,'/',['Control',13,'Drum', 
13, 'Position/Power',13,'Rate' ,13,'Positioning/Startup", 13748 
ower',13,'Profile/Heatup',13,'Rate 2',13,'(W//s)']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13, 'Rate',13, 'Positioning/Startup',13,'Power',13,'Profile/He 
atup ig, Rate 2 13." (W'/s) | 1s 

'Value','80',... 
'position',[170,140,215,160]) 





add_block('built-in/Constant',[sys,'/',['Control',13,'Drum', 
13,'Position/Power',13,'Rate',13,'Positioning/Startup',13,'P 
ower',13,'Profile/Heatup',13,'Rate 1',13,'(W//s)']]) 





set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
135 Rate 13» Positionipg/Startup' 519; "Power ' 13, Prof ile m 
atup 13 Rate 1" 13 WAS). s 

'Value','600',... 
'position',[170,265,215,285]) 





add block( '"built-in/Switch a eys, / [Control , 13 5pm 
, Position/Power',13,'Rate', 13, 'Positioning/Startup',13, 'Pow 
er',13,'Profile/Rate',13,'Mode' , 13, 'Switch']]) 

set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Ra 
Bel.13,'Mode" ,13,' Switch MIRES 

minceshold: “ss 

'position',[260,204,290,236]) 







































add block('built-in/Relay',[sys,'/',['Control',13,'Drum',13, 
'Position/Power' ,13, 'Rate',13, 'Positioning/Startup' ,13, 'Powe 
r',13,'Profile/35 kW' ,13, 'Detection']]) 

set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
E Bate',13,'Positioning/Startup',13, ‘Power’ ,13, ‘Profile/3s 
cies, Detection']],... 

'On switch value','35000',... 

'Off switch value','-35000',... 
'position',[110,208,140,232]) 





saemplock( 'built-in/Relay' ,[sys,'/',[' Control’ ,13, ‘Drum’ ,13, 
'Position/Power',13,'Rate',13,'Positioning/Startup',13,'Powe 
Eus. Profile/115 kW',13,'Detection']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
fae Hate',13,'Positioning/Startup',13, 'Power', 13, 'Profile/11 
5x 13, Detection']],... 

'On switch value','115000',... 

'Off switch value','-35000',... 

‘On output value','115000') 
set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
EuISte ,13, Positioning/Startup' , 13, ' Power ' , 13, Prof 11e/11 
DS, Detection']],... 
"position',[365,98,395,122]) 





add block('built-in/Switch',[sys, /',['Contro1',13, 'Drum' , 13 
; Position/Power',13,'Rate',13,'Positioning/Startup',13, 'Pow 
EE Profile/Gonstant',18,'Power',13, ‘Mode’ ,113, ‘Switch’ ] ] 
) 
ExwDapram([sys, /',['Control',13,'Drum', 13, 'Position/Power-, 
I hate ' ,13,' Positioning/Startup',13,'Power ', 13; "Profile/Co 
Dant 13, Power" ,13, Mode ', 13, Switch" ] 15... 

McEeshold',' 1',... 
'position',[570,204,600,236]) 
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add block ('built-in/Outport',[sys,'/',['Control',13,'Drum',1 
3,'Position/Power',13,'Rate',13,'Positioning/Startup',13,'Po 
wer',13,'Profile/Ordered',13,'Power']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13, 'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Or 
dered',13,'Power']],... 
'position',[650,210,670,230]) 


% i Subsystem 
[ Control" 13, Drum: 13, Position/Power',13,' Rate", 13, Posit 
ioning/Startup',13,'Power',13,'Profile/Order',13,'Power',13, 
'Graph']. 





new system([sys,'/',['Control',13,'Drum',193, 'Position/Power' 
,13, Rate’,13, “Positionming/Startup ,13, Power’, 13, Prof mle RE 
rder',13,'Power',13,'Graph']]) 

set _param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13, 'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Or 
der',13,'Power',13,'Graph']],'Location',[0,59,274,252]) 





add block( built-in/zInport' ,[sys, / ;['Control',135, Drums 
, Position/Power',13,'Rate',13,'Positioning/Startup',13,'Pow 
er',13,'Profile/Order',13,'Power',13,'Graph/x']]) 

set param([sys,'/*,[' Control", ,13,'Drum' 13, Position/Powe a 
13, 'Rate',13,'Positioning/Startup' ,13,/¿Power" 13, Prof 1e" 
der',13,'Power',13,'Graph/x']],... 

Positions [1655585 791) 


add_block('built-in/S-Function',[sys,'/',['Control',13,'Drum 
',13, 'Position/Power' ,13, 'Rate',13,'Positioning/Startup' ,13, 
'Power',13,'Profile/Order',13,'Power',13,'Graph/S-function', 
13, M-file which plots:,13,'lines' 13, 11) 
set _param([Sys, '/',[ ‘Control’ ,13, 'Drum*,13, ‘Position/Power 
13, 'Rate',13, 'Positioning/Startup' ,13, 'Power' ,13, 'Profile/Or 
der',13,'Power',13,'Graph/S-function',13,'M-file which 
plots 2132 mnesL l3 Jr, E 

‘function name','sfunyst') 
set param([sys,'/' [Control 513, ‘Drum’ ,13,*Position/ Powers 
13, 'Rate',13, 'Positioning/Startup',13,'Power',13,'Profillexon 
der',13,'Power',13,'Graph/S-function',13,'M-file which 
Plots 213, lines 13." 11 7e 
'Parameters','ax, color, npts, dt') 
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set param([sys, '/',['Contro1',13, 'Drum', 13, 'Position/Power', 
13,'Rate',13, 'Positioning/Startup',13, 'Power', 13, 'Profile/Or 
der',13,'Power',13,'Graph/S-function',13,'M-file which 
plots',13,/lines',13 he 

position”, (130,55, 180524315 
add line([sys,'/',['Control',13,'Drum', 13, 'Position/Power',1 
3,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Ord 
er',13,'Power',13,'Graph']],[90,65;125,65]) 
set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13,'Rate', 13, 'Positioning/Startup',13, 'Power', 13, 'Profile/Or 
der',13,'Power',13,'Graph']],... 

'Mask 


Display','plot(0,0,100,100,[83,76,63,52,42,38,28,16,11,84,11 
,11,11,90,90,11],[75,58,47,54,72,80,84,74,65,65,65,90,40,40, 
90,90])') 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13, 'Rate',13, 'Positioning/Startup',13,'Power',13, 'Profile/Or 
der',13,'Power',13,'Graph']],... 


'Mask Type','Storage scope.') 
set param([sys,'/',['Control',13,*Drum',13,'Position/Power', 
13,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Or 
der',13,'Power',13,'Graph']],... 

‘Mask  Dialogue','Storage scope using MATLAB graph 
window.\nEnter plotting ranges and line type.|Initial Time 
Range:|Initial y-min:|Initial y-max:|Storage pts.:|Line type 
fegdiow=.:xX0):' } 
set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13, 'Rate',13, 'Positioning/Startup',13, 'Power', 13, 'Profile/Or 
Een 15  Power' ,13,'Graph']],... 

‘Mask Translate','npts = @4; color = @5; ax = [0, @1, @2, 
3]; dt=-1;') 
set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13, 'Rate',13, Positioning/Startup',13, 'Power', 13, 'Profile/Or 
der',13,'Power',13,'Graph']],... 

'Mask Help','This block uses a MATLAB figure window to 
plot the input signal. The graph limits are automatically 
scaled to the min and max values of the signal stored in the 
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scope''s signal buffer. Line type must be in quotes. See the 
M-file sfunyst.m.') 
set _param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile/Or 
dendi 3 Power 13. Graph ]],... 

'Mask 
Entries','5\/-10\/10\/200\/''y-/g--/c-./w:/m*/ro/b+''\/') 


% Finished composite block 
[Control 5,13, Drum’ ‚13, Position/Power 513, Rate ',13..Posin 
ioning/Startup',13,'Power' ,13,'Profile/Order' ,13, 'Power' , 135 
'Graph']. 


set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13,'Rate',13, 'Positioning/Startup',13, 'Power', 13, 'Profile/Or 
der eis; Power’, kg, Graph |]. 

‘position’ ,[665,95,695,135] ) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Power',1 
3,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile']], 
[145,220;255,220]) 
add line([sys,'/',['Control',13,'Drum', 13, 'Position/Power',1 
3,'Rate',13,'Positioning/Startup',13, 'Power', 13, 'Profile']], 
[220,275;240,275;240,230;255,230]) 
add line([sys,'/',['Control',13, 'Drum' ,13, 'Position/Power Se 
3,'Rate',13,'Positioning/Startup',13, 'Power' , 13, 'Profile']], 
[220,150;240,150;240,210;255,210]) 
add line([sys, /',['Control',13,'Drum' ,13, 'Position/Power'j1 
3, Rate’ ,13, 'Positioning/Startup' ,13, ‘Power’ ,13, ‘Profile jie 
[355 , 2405385 ,240] ) 
add Line([sys, / ,['Contgol',13, Drumi,13, Position,Power aa 
3, 'Rate' ,13,'Positioning/Startup' ,13,'Power' ,13, ‘Profile’ ]], 
[425 ,230; 460, 230] ) 
add-line(fsys, /',['Control',13,'Drum' , 13, 'Position/Power Jom 
3j, hate ds 


'Positioning/Startup',13,'Power',13, 'Profile']],[295,220;385 
,220]) 

add line([sys,'/',['Control',13,'Drum', 13, 'Position/Power',1 
3,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile']], 
[55,20;80,20;80,220;105,220]) 
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add line([Sys, /',['Contro1',13, Drum', 13, 'Position/Power',1 
3,'Rate',13, 'Positioning/Staptup ,13, Power, 19, Profiles 
[605,220;645,220]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Power' ,1 
3.,'Rate',13, 'Positioning/Startup 13, Power ,13, Protile Vie 
[510,230;565,230]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Power',1 
3,'Rate',13,'Positioning/Startup',13,'Power',13,'Profile']], 
[400,110;540,110;540,220;565,220]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Power',1 
sEXRate',13, 'Positioning/Startup' ,13S, Power’ ,13, Profile |), 
[80,20;80,110;360,110] ) 

add line([Sys,'/',['Control',13,'Drum',13, 'Position/Power',1 
S" Bate',13,'Positioning/Startup',13, Power ', 13, 'Profile']], 
[540,210;565,210]) 

add line([Sys,'/',['Control',13,'Drum', 13, 'Position/Power',1 
3,'Rate',13,'Positioning/Startup',13, 'Power' , 13, 'Profile']], 
[618,220;640,115;660,115]) 


% Finished composite block 
['Gontrol',13,'Drum',13,'Position/Power',13,'Rate',13,'Posit 
Toning/Startup',13,'Power',13, 'Profile']. 





set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
13,'Rate',13,'Positioning/Startup',13,'Power' ,13,'Profile']] 


J 


'position',[150,91,180,144]) 














add block('built-in/Dead 
zone, |[sys,“/',[ 'Control' ,13,'Drum' ‚13, "Position/Power" , 135° 
Rate',13,'Positioning/Power Level',13,'Dead Band']]) 
set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13,'Rate',13,'Positioning/Power Level',13,'Dead Band']],.. 
"Lower value','-100',... 
supper Value','100',... 
‘position' ,[375, 47,405, 73] ) 





adsubptock't'built-in/Outport 7[sys, '/',['COnitrol' , 13, 'Drum' , 1 
3, 'Position/Power',13,'Rate',13, 'Positioning/Drum',13,'Motor 
EU 5witch']]) 

set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
Eu ute ,13,'Positioning/Drum',13,'Motor ,13, Switch']],... 
ipssition',[505,50,525,70]) 


117 


add_block('built-in/Sum',[sys,'/',['Control',13,'Drum',13,'P 
osition/Power',13,'Rate',13,'Positioning/Power',13,'Error']] 
) 
IS A IS Drums, Position/Roweri 
13,'Rate',13, ME hus. Eos 

Ln puts er 

‘position’ ‚[260, 42, 280,73]) 


add_block('built-in/Inport',[sys,'/',['Control',13,'Drum',13 
, Position/Power',13,'Rate',13,'Positioning/Rate',13,'Mode', 
13,'Selected']]) 


set param([sys,'/',['Control',13,'Drum',13, 'Position/Power', 
19,'Rate', 13, 'Positioning/Rate',13, Mode ' , 13, 'Selected']],.. 


SPORE EL C. 
'position',[40,120,60,140]) 


add block('built-in/Inport'.[sys, /',['ControL , 19, Drums 
, Position/Power',13,'Rate',13,'Positioning/Rx Power']]) 
set_param([Ssys,'/',['Control',13,'Drum',13,'Position/Power', 
13,'Rate', 13, 'Positioning/Rx Power']],... 
'position',[35,40,55,60]) 
add line([sys,'/',['Control',13,'Drum',13, 'Position/Power',1 
3, Rate "313, Positioning ]],[185,120;195,120:195,65:255, G51BI 
add_line([sys,'/',['Control',13,'Drum',13,'Position/Power',1 
3,'Rate',13,'Positioning']],[60,50;255,50]) 
add line([sys,'/',['Control',13,'Drum', 13, 'Position/Power' 
3,'Rate',13,'Positioning']],[65,130;145,130]) 
add_line([sys,'/',['Control',13,'Drum' ,13,'Position/Power' 
3,'Rate',13,'Positioning']],[410,60;500,60]) 
add_line([sys,'/',['Control',13,'Drum' ,13,'Position/Power' 
3,'Rate',13,'Positioning']],[285,60;370,60]) 
add _line([sys,'/',['Control',13,'Drum',13,'Position/Power' 
3,'Rate',13,'Positioning']],[60,50;85,50;85,105;145,105]) 
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% Finished composite block 
['Control',13,'Drum', 13; ^"Posdtyon/Power ,13, Hate IS DONE 
ionning |. 


set_param([sys,'/',['Control',13,'Drum',13,'Position/Power', 
13, 'Rate',13,'Positioning mikk. 

'position',[200,56,230,109]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position']],[510, 
200; 530,200] ) 
add_line([sys,'/',['Control',13,'Drum',13,'Position'])],[220, 
275;265,275;265,210;380,210]) 
add line([sys,'"/",[ Control’ ,13, Drum' we, ‘Position’ |], [5555 
0;95,50;95,200;110,200] ) 
add_line([sys,'/',['Control',13,'Drum',13,'Position']],[235, 
85;300,85;300,190;380,190]) 
BuSEEne([sys,'/',['Control',13,' Drum',; 13, Position ]],[1505 
200;380,200]) 
add_line([sys,'/',['Control',13,'Drum',13,'Position']],[570, 
200;605,200]) 
sadevine([sys,'/',['Control' ,13,'Drum’ ,13, Position’ ]],[645, 
200; 695,200] ) 

BgcEIne([sys, '/',['Control',13,'Drum', 13, Position ]],[645, 
200;650,360;100,360;100,275;180,275]) 


add_line([sys,'/',['Control',13, "Position’]], [655 
200;655,85;685,85]) 


dddMfne([sys,'/',['"Control',13, 'Position']],[420, 
200;470,200]) 

add_line([sys,'/',['Control',13, 'Position']],[95,5 
0;185,50;195,70]) l 
soomaime(|[sys,'/',[L'Control' ,13, 'Position']],[150, 
200;165,200;165,95;195,95]) 


% Finished composite block 
iP Control’ ,13, 'Drum',13,'Position']. 


Euwpapam([sys, /',['Control',13,'Drum', 19, "Position']1],... 
Eposrtion',[100,120,130,170]) 
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add block('built-in/To 
Workspace',[sys,'/',['Reactor',13,'Power']]) 
set param([sys,'/',['Reactor',13,'Power']],... 
'mat-name','P',... 
'buffer','40000', 
'position',[440,147,490,163]) 


% Subsystem 'Rx Power'. 


new_system([sys,'/','Rx Power' ]) 
set_param([sys,'/', ‘Rx Power'],'Location',[0,59,274,252]) 


add_block('built-in/Inport',[sys,'/','Rx Power/x']) 
set_param([sys,'/','Rx Power/x'],... 
[position | 165,55,85,, 75} ) 


add block('built-in/S-Function',[sys, / ',['Rx 
Power/S-function',13,'M-file which 
plots',13,'lines',13,'']]) 

set_param([sys,'/',['Rx Power/S-function',13,'M-file which 
plots 13, Jxnes' RE. 

'function name','sfunyst', 

'parameters','ax, color, npts, dt', 

position: [130,55 1890,75]) 
add_line([sys,'/','Rx Power'],[90,65;125,65]) 
set param([sys, /','Rx Power'],... 

' Mask 
Display','plot(0,0,100,100,[83,76,63,52,42,38,28,16,11,84,11 
,11,11,90,90,11],[75,58,47,54,72,80,84,74,65,65,65,90,40,40, 
90,90])', 

'Mask Type','Storage scope.') 

S€tlpakam( (SVs, 7 , Ax Power. |> 

‘Mask Dialogue','Storage scope using MATLAB graph 
window.\nEnter plotting ranges and line type.|Initial Time 
Range: |Initial y-min:|Initial y-max:|Storage pts.:|Line type 
(rgbw-.:xo):') 
set param([sys,'/', Rx Power’],... 

'Mask Translate','npts = @4; color = @5; ax = [0, @1, @2, 
@3]; dt=-1;') 
set param([Sys, /','Rx Power '],... 






















'Mask Help','This block uses a MATLAB figure window to plot 
the input signal. The graph limits are automatically scaled 
to the min and max values of the signal stored in the 
scope''s signal buffer. Line type must be in quotes. See the 
M-file sfunyst.m.') 

set param([sys, '/','Rx Power'],... 

Mask 

Entries','5\/-10\/10\/200\/''y-/g--/c-./w:/m*/ro/b+''\/') 


% Finished composite block 'Rx Power'. 






set_param([sys,'/','Rx Power ],... 
‘position’ ,[445,235,475,275] ) 


% Subsystem ['Temperature’',13, 'Response']. 












new system([sys,'/',['Temperature',13, 'Response']]) 
set param([sys,'/',['Temperature',13, 'Response']], 'Location' 
,[4,42,787,498]) 

open system([sys,'/',['Temperature',13,'Response']]) 


% Subsystem ['Temperature' ,13, 'Response/Tm Graph']. 





new_system([sys,'/',['Temperature' ,13, 'Response/Tm Graph' ]]) 
set_param([sys,'/',['Temperature',13,'Response/Tm 
Graph ]],'Location',[0,59,274,252]) 


add block('built-in/Inport',[sys,'/',['Temperature' , 13, ‘Resp 
onse/Tm Graph/x']]) 

set param([sys, ' /',['Temperature',13, 'Response/Tm 

px ]],... 
ipesition',[65,55,85,/5]) 

















add_block('built-in/S-Function',[sys,'/',['Temperature',13,' 
Response/Tm Graph/S-function',13,'M-file which 
poes, 13,'lines',13,'']]) 
set_param([sys,'/',['Temperature' ,13, 'Response/Tm 
Graph/S-function' ,13, 'M-file which 
sia, lines',13,'']],... 
Brumetion name’ ,'’sfunyst',... 


'parameters','ax, color, npts, dt',... 

'position',[130,55,180,75]) 
add_line([sys,'/',['Temperature',13,'Response/Tm 
Graph']],[90,65;125,65]) 
set_param([sys,'/',['Temperature',13,'Response/Tm 


Display','plot(0,0,100,100,[83,76,63,52,42,38,28,16,11,84,11 
,11,11,90,90,11],[75,58,47,54,72,80,84,74,65,65,65,90,40,40, 
9059019 4 

‘Mask Type','Storage scope. ') 
set param([sys,'/',[' Temperature',13, 'Response/Tm 
Graph']],... 


'Mask Dialogue','Storage scope using MATLAB graph 
window.\nEnter plotting ranges and line type.|Initial Time 
Range: |Initial y-min:|Initial y-max:|Storage pts.:|Line type 
(rgbw-.:x0):') 
set_param([sys,'/',['Temperature',13,'Response/Tm 
Graph']],... 

'Mask Translate','npts = @4; color = @5; ax = [0, @1, @2, 
3]; dt=-1;') 
set_param([sys,'/',['Temperature' ,13, 'Response/Tm 
Graph']],... 

'Mask Help','This block uses a MATLAB figure window to 
plot the input signal. The graph limits are automatically 
scaled to the min and max values of the signal stored in the 
scope''s signal buffer. Line type must be in quotes. See the 
M-file sfunyst.m.') 
set_param([sys,'/',['Temperature' ,13, 'Response/Tm 
Graph']],... 

'Mask 
Entries','5\/-10\/10\/200\/''y-/g--/c-./w:/m*/ro/bt''\/') 


% Finished composite block ['Temperature',13,'Response/Tm 
Graph']. 


set param([sys, ' /',['Temperature',13, 'Response/Tm 
Graph jl,- 
'position',[625,315,655,355]) 





add block('built-in/Note',[sys, '/',[' Temperature', 13, 'Respon 

se/Tm']]) 

set_param([sys,'/',['Temperature',13,'Response/Tm']],... 
'position',[300,255,305,260]) 


% Subsystem ['Temperature',13, 'Response/Tm Dynamics']. 


new system([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics']]) 

set param([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics']],'Location',[18,61,780,513]) 


add block('built-in/Mux',[sys, /',['Temperature',13, 'Respons 
e/Tm Dynamics/Tm Op']]) 
set param([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics/Tm Op']],... 

EUiputs','2',... 

'position',[300,231,330,264]) 


add block('built-in/Switch',[sys,'/',['Temperature' , 13, 'Resp 
onse/Tm Dynamics/TFE',13,'Loading']]) 
set param([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics/TFE',13,'Loading']],... 

‘Threshold','30000',... 

gposs$tion',{220,239,250,271]) 


add block('built-in/Outport',[sys, /',['Temperature' , 13, 'Res 
ponse/Tm Dynamics/Tm Rate']]) 
set param([sys,'/',['Temperature' ,13, 'Response/Tm 
Dynamics/Tm Rate']],... 

gsSsrtion',[695,287,715,313]) 


add block('built-in/Switch',[sys,'/',['Temperature' , 13, 'Resp 
onse/Tm Dynamics/Cs add',13,' TFE Load']]) 
set param([sys,'/',['Temperature',13, 'Response/Tm 
Pamics/Cs add',13,'TFE Load']],... 

ampeshold','Teos',... 

'position',[605,284,635,316]) 
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add block('built-in/Inport',[sys, '/',['Temperature' ,13, 'Resp 
onse/Tm Dynamics/Time']]) 

set param([sys,'/',['Temperature' , 13, 'Response/Tm 
Dynamics/Time']],... 
PROB 2 x. 

position", ,[505,290,525,310]) 




































add_block('built-in/MATLAB 
Fcn',[sys,'/',['Temperature' ,13,'Response/Tm Dynamics/Tm 
Equilibrium',13,'(Low Power)']]) 

set param([sys,'/',['Temperature' , 13, 'Response/Tm 
Dynamics/Tm Equilibrium',13,'(Low Power)']],... 

'MATLAB Fcn','298.09755*0.012550343*u' ,... 

"Output Width we a 
position ;[1105314,185,356]) 






add block('built-in/MATLAB 
Fen',[sys,'/',[' Temperature’ ,13, 'Response/Tm Dynamics/Tm 
Equilibrium',13,'(High Power)']]) 
set_param([sys,'/',['Temperature',13,'Response/Tm 
Dynamics/Tm Equilibrium',13,' (High Power)']],... 

‘MATLAB Fen' ,'664.13298+1 .238966e-8*u*u',... 

Output Width ee: 
'position',[100,119,175,161]) 






add block('built-in/MATLAB 
Fen',[sys,'/',['Temperature' ,13, 'Response/Tm 
Dynamics/Startup Dynamics']]) 
set param([sys,'/',['Temperature' ,13, 'Response/Tm 
Dynamics/Startup Dynamics']],... 

' MATLAB Fen','(6.3617662 -0.0056230422*u(1)/u(2)*825 
:11379988541*u(2)70(1)7/825)*u(2)7825 50€ 
Output Width Ti. 
'position',[390,334,465,376]) 






add block('built-in/MATLAB 
Fcn',[sys,'/',['Temperature',13,'Response/Tm 
Dynamics/Operating',13,'Dynamics']]) 
set param([SyS, '/',['Temperature',13, 'Response/Tm 
Dynamics/Operating',13,'Dynamics']],... 
IMATIEAB (CN, 0002680 -30(2)8) T, . 2 
poueputoWasdtb 1... 
'position',[390,229,465,271]) 
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add block('built-in/Inport',[sys, /',[' Temperature ',13, 'Resp 
onse/Tm Dynamics/Power']]) 
set param([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics/Power']],... 

Port. <3' se 

'position',[20,130,40,150]) 


add_block('built-in/Inport',[sys,'/',['Temperature',13, Resp 
onse/Tm Dynamics/Tm']]) 
set param([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics/Tm']],... 

Eorientation',2,... 

'position',[590,205,610,225]) 
add line([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics']],[335,250;385,250]) 
add line([sys,'/',['Temperature',13, 'Response/Tm 
Dynamics']],[585,215;275,215;275,240;295,240]) 
add_line([sys,'/',['Temperature',13,'Response/Tm 
Dynamics']],[255,255;295,255]) 
add line([sys,' /',['Temperature' , 13, 'Response/Tm 
Dynamics']],[180,140;205,140;215,245]) 
add line([sys,'/',['Temperature' , 13, 'Response/Tm 
Dynamics']],[530,300;600,300]) 
add_line([sys,'/',['Temperature',13,'Response/Tm 
Dynamics']],[640,300;690,300]) 
add line([sys,'/',['Temperature' , 13, 'Response/Tm 
Dynamics']],[470,355;565,355;565,310;600,310]) 
add_line([sys,'/',['Temperature',13,'Response/Tm 
Dynamics']],[45,140;60,140;60,255;215,255]) 
add line([sys,'/',['Temperature' ,13, 'Response/Tm 
Dynamics']],[190,335;190,265;215,265] ) 
add line([sys,'/',['Temperature' ,13, 'Response/Tm 
Dynamics']],[60,145;60,140;95, 140] ) 
add line([sys, /',['Temperature', 13, 'Response/Tm 
Dynamics']],[60,255;60,335:105,335]) 
add _line([sys,'/',['Temperature',13,'Response/Tm 
Dynamics']],[470,250;565,250;565,290;600,290]) 
add_line([sys,'/',['Temperature',13,'Response/Tm 
eynamues']),[335,250;355, 250; 355, 355;385,355] ) 





% Finished composite block ['Temperature',13,'Response/Tm 
Dynamics']. 


set param([sys, '/',['Temperature',13, 'Response/Tm 
Dynamics']],... 
'position',[300,307,330,363]) 


add _block('built-in/Integrator',[sys,'/',['Temperature',13,' 
Response/Tm Integrator']]) 
set param([sys,'/',['Temperature',13, 'Response/Tm 
Integrator']],..- i 
"Initial TO ser 
'position',[410,322,440,348]) 


add_block('built-in/Note',[sys,'/',['Temperature' ,13, 'Respon 

se/Tf']]) 

set param([sys,'/',['Temperature',13, 'Response/Tf ']],... 
'position',[290,105,295,110]) 


add_block('built-in/To 
Workspace',[Sys,'/',['Temperature',13,'Response/time']]) 


set param([sys,'/',['Temperature',13, 'Response/time']],... 
'mat-name','t',... 
'puffer','40000',... 
'position',[210,227,260,243]) 


add block('built-in/Integrator',[sys,'/',['Temperature' ,13,' 
Response/Tf Integrator']]) 
set param([sys,'/',['Temperature',13, 'Response/Tf 
iIntegrator il, : 

ANA TLAL Ton. 

'position',[400,162,430,188]) 


% Subsystem ['Temperature',13,'Response/Tf Dynamics']. 


new system([sys,'/',['Temperature', 13, 'Response/Tf 
Dynamics']]) 
set_param([sys,'/',[' Temperature’ ,13, Response/Tf 
Dynamics']],'Location',[118,104,709,472]) 





add block( 'built-in/OutporpGs (sys, / ,{ Temperature T3 RES 
ponse/Tf Dynamics/Tf Rate']]) 

set param([Sys,'/',['Temperature' , 13, 'Response/Tf 
Dynamics/Tf Rate']],... 
'position',[545,150,565,170]) 













add, block('built-in/Sum',[sys,'/',['Temperature' , 13, 'Respons 
e/Tf Dynamics/Sum']]) 

set param([sys,'/',['Temperature' ,13, 'Response/Tf 
Dynamics/Sum']],... 
'position',[410,285,430,305]) 


add block('built-in/Switch',[sys,'/',['Temperature',13,'Resp 
onse/Tf Dynamics/TFE',13,'Load',13,'Switch']]) 

set _param([sys,'/',['Temperature',13,'Response/Tf 
bimiamies/TFE',13,'Load',13,'Switch']],... 
'Threshold','Tload',... 
'position' ,[485,144,515,176] ) 






















% Subsystem ['Temperature' ,13, 'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']. 





new system([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Cs 
Addition',13,'Dynamics']],'Location',[36,179,729,468]) 


add block('built-in/Switch',[sys, '/',['Temperature', 13, 'Resp 

onse/Tf Dynamics/Cs Addition',13,'Dynamics/Switch']]) 

set param([sys,'/',['Temperature' , 13, 'Response/Tf 

Dynamics/Cs Addition',13,'Dynamics/Switch']],... 
gnreshold','Tcs',... 

'position',[160,29,190,61]) 











add block('built-in/Inport',[sys, '/',['Temperature' ,13, 'Resp 
onse/Tf Dynamics/Cs Addition',13,'Dynamics/Time']]) 
set param([Sys, ' /',['Temperature',13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Time']],... 

ort, '2',... 
geos#tion',[15,5,395725]) 
















add block('built-in/Switch',[sys, '/',['Temperature' , 13, 'Resp 
onse/Tf Dynamics/Cs Addition',13,'Dynamics/Cs 
Addition',13,'Switch']]) 

set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Cs 
Addition',13,'Switch']],... 
'Threshold','Tes',... 
'position',[585,94,620,136]) 





add block('built-in/Outport',[sys, '/', ['Temperature' , 13, 'Res 
ponse/Tf Dynamics/Cs Addition',13,'Dynamics/Cs 
Heatup',13,'Rate']]) 

set _param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Cs 
Heatup',13, 'Rate']],... 
'orientation',2,... 
'position',[580,15,600,35]) 























add_block('built-in/Constant',[sys,'/',['Temperature',13,'Re 
sponse/Tf Dynamics/Cs Addition',13,'Dynamics/Constant']]) 
set _param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Constant']],... 
'orientation',2,... 

A 
'position',[600,185,620,205]) 






add block('built-in/MATLAB 
Fon',[sys,'/',['Temperature' ,13, 'Response/Tf Dynamics/Cs 
Addition' ,13, 'Dynamics/Tf Dynamics' ,13,'Cs Addition']]) 
set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Tf Dynamics',13,'Cs 
Acton n 

' MATLAB Fen','sqrt(polyval([-7.6606914e-6 0.032593747 
0.0179884],u(1)))*u(2)/115e3',... 

‘Output Width’ , ww) 
set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Tf Dynamics',13,'Cs 
Addition']],... 
'position',[445,79,520,121]) 


add _block('built-in/Memory',[sys,'/',['Temperature',13,'Resp 
onse/Tf Dynamics/Cs Addition',13,'Dynamics/Memory']]) 

set param([sys, '/',['Temperature' ,13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Memory']],... 


'position',[220,30,260,60]) 


add _block('built-in/Mux',[sys,'/',['Temperature' ,13,'Respons 
e/Tf Dynamics/Cs Addition',13,'Dynamics/Mux']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Mux']],... 

Zinputs' 727, 


'position',[395,61,430,139]) 


add block('built-in/Abs',[sys,'/',[' Temperature' , 13, 'Respons 

e/Tf Dynamics/Cs Addition',13,'Dynamics/Abs']]) 

set param([sys,'/',['Temperature',13, 'Response/Tf 

Dynamics/Cs Addition',13,'Dynamics/Abs']],... 
'position',[350,68,380,92]) 


add block('built-in/Sum',[sys, '/',['Temperature' , 13, 'Respons 
e/Tf Dynamics/Cs Addition',13,'Dynamics/Sum']]) 
set param([sys,'/',['Temperature',13, 'Response/Tf 
Pees Cs Addition' ,13, 'Dynamics/Sum']],... 
Enmbuts ga-t', 
‘position’ ‚[290, ‚64, 320,91]) 


add block('built-in/Demux',[Sys,'/',['Temperature', 13, 'Respo 
nse/Tf Dynamics/Cs Addition',13,'Dynamics/Demux']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Cs Addition' ,13, 'Dynamics/Demux']],... 

EUEDUtS' ,'2',... 

ipngsrtion' ,[50,75,105,110]) 


add block('built-in/Note',[sys, /',['Temperature' , 13, 'Respon 

se/Tf Dynamics/Cs Addition',13,'Dynamics/Tf']]) 

set param([Sys,'/',['Temperature', 13, 'Response/Tf 

Dynamics/Cs Addition',13,'Dynamics/Tf']],... 
'position',[75,45,80,50]) 


add_block('built-in/Note',[sys,'/',['Temperature',13,'Respon 

se/Tf Dynamics/Cs Addition',13,'Dynamics/P']]) 

set _param([sys,'/',['Temperature',13,'Response/Tf 

Dynamics/Cs Addition',13,'Dynamics/P']],... 
Eposstuion',[130,105,135,110]) 


add block('built-in/Inport',[sys, '/',[' Temperature', 13, ' Resp 
onse/Tf Dynamics/Cs Addition',13,'Dynamics/State']]) 








set param([sys,'/',['Temperature' , 13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/State']],. 
'position',[10,85,30,105]) 




































add block('built-in/Note',[sys, /',['Temperature' , 13, 'Respon 
se/Tf Dynamics/Cs Addition',13,'Dynamics/Tf initial']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics/Tf initial']],... 
Position", [310,30,315,35]/) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[525,100;580,100]) ' 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[595,195;560,195;560,130;580,130]) 





add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[625,115;635,115;635,25;605,25]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Aadalition” , 13, Dynamics”]],[40, 153550015 9390071535803 1151) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[195,45;215,45]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[265,45;265,25;140,25;140,35;155,3 
5]) 

add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[35,95;45,95]) 

add line([sys,'/',['Temperature',13, 'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[435,100;440,100]) 

add line([Sys,'/',['Temperature',13, 'Response/Tf Dynamics/Cs 
Addition , 13, Dynamics’ |i NO 8510085160 55;155 5511 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[95,15;95,45;155,45]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[110,100;245,100;245,120;390,120]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[265,45;275,45;285,70]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[110,85;285,85]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[385,80;390,80]) 
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add line([sys,'/',['Temperature', 13, 'Response/Tf Dynamics/Cs 
Addition',13,'Dynamics']],[325,80;345,80]) 










% Finished composite block ['Temperature',13,'Response/Tf 
Dynamics/Cs Addition' ,13,'Dynamics']. 

















set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Cs Addition',13,'Dynamics']],... 
'position',[255,187,285,238]) 






add block('built-in/MATLAB 
Fen',[sys,'/',['Temperature',13,'Response/Tf Dynamics/Tf 
Dynamics',13,'(No Loading)']]) 
set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Tf Dynamics',13,' (No Loading)']],... 

' MATLAB Fcn','1.201131601e-4*u(2) *polyval([-2.0220458e-8 
7.3516563e-5 -0.087210261 20.092539996],u(1))',... 

{output Width','1') 
set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Tf Dynamics',13,'(No Loading)']],... 
'position',[195,279,270,321]) 


add block('built-in/Inport',[sys, '/',['Temperature' , 13, 'Resp 
onse/Tf Dynamics/Time']]) 

set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Time']],... 

Seem ,'3S',... 
gooSttion',[15}150,35,170]) 

















% Subsystem ['Temperature',13,'Response/Tf | Dynamics/Tf 
Dynamics',13,'(Operating)']. 





new_system([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf 
Dynamics',13,'(Operating)']],'Location',[62,71,721,547]) 


add_block('built-in/Inport',[sys,'/',['Temperature',13,'Resp 
onse/Tf Dynamics/Tf Dynamics',13,'(Operating)/Power']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Power']],... 
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uo 22 Sects 
'position',[35,175,55,195]) 


add block('built-in/MATLAB 
Fcn',[sys, '/',['Temperature' , 13, 'Response/Tf Dynamics/Tf 
Dynamics',13,'(Operating)/Tf Equilibrium']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Tf Equilibrium']],... 
'MATLAB Fcn','polyval( [ 2.7674712e-13 O -0.019137039 
299.54525] , U ) + 9.9191989*sqrt(u)', 
"Output Width 1) : 
set param([sys,'/',[' Temperature ,13, 'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Tf Equilibrium']],... 
'position',[105,164,180,206]) 


add block('built-in/Sum',[sys,'/',['Temperature', 13, 'Respons 
e/Tf Dynamics/Tf 
Dynamics',13,'(Operating)/Tf',13,'Deviance']]) 
set param([sys,'/',['Temperature', 13, 'Response/Tf 
Dynamics/Tf 
Dynamics”, SS agas s EAE 13,'Deviance']],... 

'inputs', i 


‘position’ ‚[225, 114, 250,141]) 


add block('built-in/Constant',[sys,'/',['Temperature' , 13, ' Re 
sponse/Tf Dynamics/Tf Dynamics' ,13, ' (Operating) /Heatup 
Rate']]) 
set_param([sys,'/',['Temperature' ,13, 'Response/Tf 
Dynamics/Tf Dynamics' ,13,'(Operating) /Heatup Rate']],... 
'Value','-0.022569155',... 
'position' ,[215,302,365, 328] ) 


add_block('built-in/Constant',[sys,'/',['Temperature',13,'Re 
sponse/Tf  Dynamics/Tf  Dynamics',13,'(Operating)/Cooldown 
Rate']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Cooldown Rate']],... 
'Value','-0.01',.. 
‘position’ ,[295,191,435,219] ) 


add block('built-in/Switch',[sys,'/',['Temperature ' , 13, 'Resp 
onse/Tf Dynamics/Tf Dynamics',13,' (Operating)/Switch']]) 






























set param([sys, '/',['Temperature',13, 'Response/Tf 
Dynamics/Tf Dynamics' ,13,' (Operating) /Switch']],... 
‘position' ,[460,234,490, 266] ) 


add_block('built-in/Product',[sys,'/',['Temperature',13,'Res 
ponse/Tf Dynamics/Tf Dynamics',13,'(Operating)/Tf 
rate',13,'Calculation']]) 

set param([sys,'/',['Temperature', 13, 'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Tf 
fame ,13,' Calculation) ],..- 
position! ,[525,123,555,147]) 


add _block('built-in/Outport',[sys,'/',['Temperature',13,'Res 
ponse/Tf Dynamics/Tf Dynamics',13,'(Operating)/Tf Rate']]) 
set _param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Tf Rate']],... 
8position',[615,125,635,145]) 


adggeblock('built-in/Inport',[sys,' / ',['Temperature ,13, Resp 
onse/Tf Dynamics/Tf Dynamics',13,'(Operating)/Tf']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)/Tf']],... 
'position',[30,110,50,130]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Tf 
Dynamics',13,'(Operating)']],[60,185;100,185]) 
add_line([sys,'/',['Temperature',13, 'Response/Tf Dynamics/Tf 
Dynamics',13,'(Operating)']],[185,185;195,185;195,135;220,13 
5]) 
add line([sys,'/',['Temperature', 13, ' Response/Tf Dynamics/Tf 
Dynamics' ,13,' (Operating) ']],[55,120;220,120] ) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Tf 
Dynamics' ,13,' (Operating) ']],[495,250;500 ,250;500, 140;520,14 
0]) 
add line([sys,'/',['Temperature' , 13, 'Response/Tf Dynamics/Tf 
Dynamics' ,13,' (Operating) ']],[560,135;610,135]) 
add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Tf 
Dynamics' ,13,' (Operating) ']],[255,130;275,130;275,250; 455,25 
0]) 
add line([sys,'/',['Temperature', 13, 'Response/Tf Dynamics/Tf 
Dynamics',13,'(Operating)']],[440,205;445,205;455,240]) 
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add line([sys,'/',['Temperature',13, 'Response/Tf Dynamics/Tf 
Dynamics',13,'(Operating)']],[370,315;390,315;390,260;455,26 
01) 

add_line([sys,'/',['Temperature',13,'Response/Tf Dynamics/Tf 
Dynamics',13,'(Operating)']],[275,130;520,130]) 










% Finished composite block 
Dynamics/Tf Dynamics',13, 


['Temperature',13, 'Response/Tf 
' (Operating) ']. 














set param([Sys,'/',['Temperature',13, 'Response/Tf 
Dynamics/Tf Dynamics',13,'(Operating)']],... 
"position ,[p295,50,325, 103] 


add block('built-in/Inport',[sys, '/',['Temperature' , 13, 'Resp 
onse/Tf Dynamics/Tf']]) 















set param([sys, /',['Temperature',13, 'Response/Tf 
Dynamics/Tf']],... 
'position',[20,20,40,40]) 





add_block('built-in/Mux',[sys,'/',['Temperature',13,'Respons 
e/Tf Dynamics/Mux']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf 
Dynamics/Mux']],... 
NOUS Le er 
‘position , (105,101,735, 134) 


add block('built-in/Inport',[sys,'/',[' Temperature ',13, 'Resp 
onse/Tf Dynamics/Power']]) 
set param([sys, /',['Temperature',13, 'Response/Tf 
Dynamics/Power']],... 

A see 

positions 1M5,80,35/1001!) 
add_line([sys,'/',['Temperature',13,'Response/Tf 
Dynamics']],[140,120;160,120;160,300;190,300]) 
add line([Sys,'/',['Temperature',13, 'Response/Tf 
Dynamics']],[160,200;250,200]) 
add line([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics']],[435,295;450,295;450,170;480,170]) 
add_line([sys,'/',['Temperature',13,'Response/Tf 
Dynamics']],[290,215;375,215;375,290;405,290]) 
















134 





add line([sys,'/',['Temperature ,13, 'Response/Tf 
Dynamics']],[275,300;405,300]) 

add line([Sys,'/',['Temperature',13, 'Response/Tf 
Dynamics']],[330,80;370,80;370,150;480,150]) 

add line([sys,'/',['Temperature' ,13, 'Response/Tf 
Dynamics']],[45,30;65,30;65,110;100,110]) 
add_line([sys,'/',['Temperature',13,'Response/Tf 
Dynamics']],[40,160;480,160]) 

add line([Sys,'/',['Temperature', 13, 'Response/Tf 
Dynamics']],[520,160;540,160]) 
add_line([sys,'/',['Temperature',13,'Response/Tf 
Dupamics']],[75,160:75,225:250,225]) 

add line([sys,'/',['Temperature' , 13, 'Response/Tf 
Dynamics']],[65,65;290,65]) 
add_line([sys,'/',['Temperature',13,'Response/Tf 
Dynamics']],[40,90;290,90]) 
add_line([sys,'/',['Temperature',13,'Response/Tf 
Dynamics']],[50,90;50,125;100,125]) 































% Finished composite block ['Temperature',13, 'Response/Tf 
Dynamics']. 





set param([sys,'/',['Temperature',13, 'Response/Tf 
Dynamics']],... 
'position',[300,148,330,202]) 











add block('built-in/To 
Workspace',[sys,'/',['Temperature' ,13, 'Response/Tm out']]) 


set_param([sys,'/',['Temperature',13,'Response/Tm out']],... 
'mat-name','Tm',... 
‘buffer','40000',... 
‘position’ ,[680, 392,730,408] ) 





add_block('built-in/Mux',[sys,'/',['Temperature',13,'Respons 
e/Mux']]) 
set_param([sys,'/',['Temperature' ,13, 'Response/Mux']],... 
ENDBUtS','3',... 
'position',[685,244,715,276]) 


add block('built-in/MATLAB 
Fen',[Sys,'/',['Temperature',13, 'Response/Te Calculation']]) 
set param([sys,'/',['Temperature',13, 'Response/Te 
Calculation']],.. 
"MATLAB Fcn','-6.8425844*0.99243093*u' ,... 
"position',[555,181,630,219]) 


% Subsystem ['Temperature',13,'Response/ Tf Graph']. 


new system([sys, '/',['Temperature' ,13, 'Response/ Tf 
Graph']]) 
set param([sys,'/',['Temperature' , 13, 'Response/ Tf 


Graph']],'Location',[0,59,274,252]) 


add block('built-in/Inport',[sys, '/',['Temperature ' , 13, 'Resp 
onse/ Tf Graph/x']]) 
set param([sys, /',['Temperature',13, 'Response/ Tf 
Grape ll... 

POSItion' 165,55, 85 751) 


add_block('built-in/S-Function',[sys,'/',['Temperature',13,' 


Response/ Tf Graph/S-function',13,'M-file which 
DICES 13: lines 13553010 

set param([sys,'/',['Temperature',13, 'Response/ Toh 
Graph/S-function',13,'M-file which 


plots 137 lines 13.' qs .. 
‘function name','sfunyst',... 
'parameters','ax, color, npts, dt',... 
'position',[130,55,180,75]) 


add_line([sys,'/',['Temperature',13,'Response/ Tf 
Graph']],[90,65;125,65]) 
set_param([sys,'/',['Temperature’' ,13, 'Response/ Ti 
Graph']],... 

'Mask 


Display','plot(0,0,100,100,[83,76,63,52,42,38,28,16,11,84,11 
,11,11,90,90,11],[75,58,47,54,72,80,84,74,65,65,65,90,40,40, 
99.90 ey... 

‘Mask Type','Storage scope.') 
set _param([sys,'/',['Temperature' ,13,'Response/ Tf 
Gap ve 

‘Mask Dialogue','Storage scope using MATLAB graph 
window.\nEnter plotting ranges and line type.|Initial Time 
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Range:|Initial y-min:|Initial y-max:|Storage pts.:|Line type 
(rgbw-.:xo): ') 
set param([Sys, '/',['Temperature',13, 'Response/ If 
Graph']],... 

'Mask Translate','npts = @4; color = @5; ax = [0, @1, @2, 
@3]; dt=-1;') 


set_param([sys,'/',['Temperature',13,'Response/ Tf 
Gkaph']|],... 

‘Mask Help','This block uses a MATLAB figure window to 
plot the input signal. The graph limits are automatically 
scaled to the min and max values of the signal stored in the 
scope''s signal buffer. Line type must be in quotes. See the 
M-file sfunyst.m.') 
set param([sys,'/',['Temperature', 13, 'Response/ Tf 
asaph 1],... 

Mask 
Enbp3es','5V/-10V/101/200V/' ' y-/g- -/C- ./w:/m*/ro/bt' 'NX/') 


$ Finished composite block ['Temperature',13,'Response/ Tf 
Graph']. 


set param([sys,'/',['Temperature',13, 'Response/ ET 
Graph']],... 
'position',[475,210,505,250]) 


add block('built-in/To 
Workspace',[sys,'/',['Temperature',13,'Response/Tf Out']]) 
set_param([sys,'/',['Temperature',13,'Response/Tf Out']],... 

‘mat-name’,'Tf',... 

"buffer','40000',... 

'position',[680,167,730,183]) 


add block('built-in/Clock',[sys,'/',[' Temperature ' , 13, 'Respo 

nse/Clock']]) 

set param([Ssys, /',['Temperature',13, 'Response/Clock']],... 
'position',[140,185,160,205]) 


add block('built-in/Inport',[sys, /',[' Temperature ', 13, 'Resp 
onse/Power']]) 


dior 


set param([sys,'/',['Temperature',13,'Response/Power']],... 


Position, [35,165,55,185]) 


add block('built-in/Outport',[sys,'/ 
ponse/Trx']]) 


',[ Temperature ',13, 'Res 


set param([sys,'/',['Temperature',13, 'Response/Trx']],... 


ipesityon »[750:250, 770,270] ) 

add_line([sys,'/',['Temperature',13, 
5,195]) 

add line([sys, /',['Temperature',13, 
5.335: 295335] 
add_line([sys,'/',['Temperature',13, 
529351) 

add line([sys,'/',['Temperature',13, 


'Response']1],[165,195;29 


"Response 1], [175,195; 17 


'Response']],[175,235;20 


'Response']],[435,175;44 


5,175;445,135;280,135;280,155;295,155]) 


add line([Sys, '/',['Temperature',13, 
0,175;440,230;470,230]) 
add_line([sys,'/',['Temperature',13, 
5,285;280,285;280,315;295,315]) 
add_line([sys,'/',['Temperature',13, 
0,335]) 
add_line([sys,'/',['Temperature',13, 
0,400;675,400]) 

add line([Sys, '/',['Temperature' ,13, 
27530 | } 
add_line([sys,'/',['Temperature' ,13, 
55175) 

add line([SsyS, '/',['Temperature' ,13, 
,175]) 


add line([Sys,'/',['Temperature' ,13, 
5,260]) 

add line([sys,'/',['Temperature',13, 
5,270;680,270]) 
add_line([sys,'/',['Temperature' ,13, 
0,200;650,260;680, 260] ) 
add_line([sys,'/',['Temperature' ,13, 
5.175.665, 250:6807,250]) 

add line([sys,'/',['Temperature' ,13, 
5,200;550,200]) 

add line([sys, /' ,['Temperature 139, 
5.172510 
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'Response']],[435,175;44 
‘Response ']],[445,335;44 
'Response']],[445,335;62 
'Response']],[590,335;59 
'Response']],[335,335;40 
'Response']],[335,175;39 


‘Response ']],[60,175;295 


'Response']],[720,260;74 
'Response']], [545,335;54 
'Response']],[635,200;65 
'Response']],[435,175;66 
'Response']],[525,175;52 


'Response']],[660,175;67 





add line([sys, '/',['Temperature',13, 'Response']],[90,175,90, 
355;295,355]) 


% Finished composite block ['Temperature',13,'Response']. 


set _param([sys,'/',['Temperature',13,'Response']],... 
'position',[540,70,570,120]) 


add block('built-in/Note',[sys,'/',['Temperature', 13, 'Feedba 

ck']]) 

set _param([sys,'/',['Temperature',13,'Feedback']],... 
position [15051551019 


% Subsystem ['Reactivity',13,'Model']. 


new_system([sys,'/',['Reactivity',13,'Model']]) 
set_param([sys,'/',['Reactivity',13,'Model']],'Location',[20 
210/5687 ,518]) 


add_block('built-in/Inport',[sys,'/',['Reactivity',13,'Model 

PT Rx”) )) : 

set param([sys,'/',['Reactivity',13, 'Mode1/T Rx']],... 
iposition',[45,85,65,105]) 


add block('built-in/Inport',[sys,'/',['Reactivity',13,'Model 
(Ogun ,13,'Position”]]) 
set param([sys,'/',['Reactivity',13, Model/Drum',13, 'Positio 
mob... 

Et, 2',..— 

'position',[40,175,60,195]) 


add_block('built-in/MATLAB 
Fen',[sys,'/',['Reactivity',13,'Model/Drum Reactivity']]) 
set _param([sys,'/',['Reactivity',13,'Model/Drum 
Heactivity']],... 

'MATLAB Fcn','polyval([6.89e-13 -2.33e-10 3.28e-9 4.57e-6 
-5.88e-5 -0.034009],u)',... 

output Width'.'1',... 

‘position’ ,[205,166 ,290,204] ) 
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% Subsystem 
['Reactivity' ,13, 'Model/Temperature' ,13, 'Reactivity']. 


new_system([sys,'/',['Reactivity' ,13, 'Model/Temperature' ,13, 
'Reactivity']]) 

set param([sys, ' /',['Reactivity',13, Model/Temperature',13,' 
Reactivity']],'Location',[45,158,708,527]) 


add _block('built-in/Demux',[sys,'/',['Reactivity',13,'Model/ 
Temperature',13,'Reactivity/Demux']]) 
set param([sys,'/',['Reactivity',13, 'Model/Temperature',13,' 
Reactivity/Demux']],... 

OUTPUTS ge 

'position',[110,154,165,186]) 


add block('built-in/Inport',[sys, '/',['Reactivity',13,' Model 
/Temperature',13, 'Reactivity/Reactor',13,'Temperatures']]) 
set param([sys,'/',['Reactivity',13, 'Model/Temperature',13,' 
Reactivity/Reactor',13,'Temperatures']],... 


'position',[50,160,70,180]) 


add_block('built-in/MATLAB 
Fen',[Sys,'/',['Reactivity',13,'Model/Temperature',13,'React 
ivity/Fuel',13,'Reactivity (Tf)']]) 
set_param([sys,'/',['Reactivity',13,'Model/Temperature',13,' 
Reactivity/Fuel',13,'Reactivity (Tf)']],... 
'MATLAB Fcn','(0.432e-5* (u- To) -6.4e-3* (sqrt(u/To)-1)) ' ,... 
'position',[230,246,305,284]) 


add block('built-in/MATLAB 
Fcn',[sys,'/',['Reactivity',13,'Model/Temperature',13,'React 
ivity/Emitter',13,'Reactivity (Te)']]) 

set param([sys,'/',['Reactivity',13,'Model/Temperature' ,13, ' 
Reactivity/Emitter',13,'Reactivity (Te)']],... 

' MATLAB Fcn','(8.52e-2-4.26e-2*( 
4.0257624+0.026717601*sqrt (u) -44.297412/sqrt(u) 
))*.01-4.93411008E-05') 
set param([sys,'/',['Reactivity',13,'Model/Temperature',13,' 
Reactivity/Emitter',13,'Reactivity (Te)']],... 

'position',[225,151,300,189])) 
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add block('built-in/MATLAB 
Fcn',[sys, '/',['Reactivity',13, ' Model/Temperature',13, 'React 
ivity/Moderator',13,'Reactivity (Tm)']]) 
set param([sys,'/',['Reactivity',13,'Model/Temperature',13,' 
Reactivity/Moderator',13,'Reactivity (Tm)']],... 

'MATLAB Fcn','polyval([-8.22e-14 1.6e-10 -1.11e-7 2.92e-5 
1576-3017. U:T0)* 401 cas 

'position',[225,56,300,94]) 


add block('built-in/Sum',[sys, '/',['Reactivity',13, Model/Te 
mperature',13,'Reactivity/Sum']]) 


set param([sys,'/',['Reactivity',13, Model/Temperature',13, 
Reactivity/Sum']],... 

BEDDUtS', ttt',... 

'position',[415,152,435,188]) 


add block('built-in/Outport',[sys, '/',['Reactivity',13, 'Mode 
1/Temperature',13,'Reactivity/Temperature',13,'Reactivity']] 
) 
set _param([sys,'/',['Reactivity',13,'Model/Temperature',13,' 
Reactivity/Temperature',13,'Reactivity']],... 
'position',[535,160,555,180]) 
add_line([sys,'/',['Reactivity',13,'Model/Temperature',13,'R 
eactivity']],[440,170;530,170]) 
add_line([sys,'/',['Reactivity',13,'Model/Temperature',13,'R 
eẹactivity']],[305,75;355,75;355,160;410,160)) 
add_line([sys,'/',['Reactivity',13,'Model/Temperature',13,'R 
eactivity']],[305,170;410,170]) 
add _line([sys,'/',['Reactivity',13, 'Model/Temperature',13,'R 
eactivity']],[75,170;105,170]) 
add_line([sys,'/',['Reactivity',13, 'Model/Temperature',13,'R 
pactryity']],[170,170;220,170]) 
add_line([sys,'/',['Reactivity',13,'Model/Temperature',13,'R 
eactivity']],[170,160;180,160;180,75;220,75]) 
add_line([sys,'/',['Reactivity',13,'Model/Temperature',13,'R 
eactivity']],[170,180;180,180;180,265;225,265]) 
add_line([sys,'/',['Reactivity',13,'Model/Temperature',13,'R 
eactivity']],[310,265;355,265;355,180;410,180]) 
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% Finished composite block 
['Reactivity',13, 'Model/Temperature' ,13, 'Reactivity']. 


set_param([{sys,'/',['Reactivity' ,13, 'Model/Temperature’ ,13, ' 
Reactivity']],... 
'position',[230,69,260,121]) 


add_block('built-in/Note',[sys,'/',['Reactivity',13,'Model/R 
x',13,'Temperatures']]) 
set_param([sys,'/',['Reactivity',13,'Model/Rx',13,'Temperatu 
res Je 

‘position: [1655071710 551) 


add block('built-in/Sum',[sys,'/',['Reactivity',13, Mode1/Su 

m']]) 

set_param([sys,'/',['Reactivity',13,'Model/Sum']],... 
'position',[365,125,385,145]) 


add_block('built-in/Note',[sys,'/',['Reactivity' ,13, 'Model/D 

rum Angle']]) 

set param([sys,'/',['Reactivity',13,'Model/Drum Angle']],... 
'position',[140,140,145,145]) 


add_block('built-in/To 
Workspace',[sys,'/',['Reactivity' ,13, 'Model/Temp' ,13, 'Reacti 
VO 
set_param([sys,'/',['Reactivity',13,'Model/Temp',13,'Reactiv 
IV 

'mat-name','r T',... 

'buffer','40000',... 

'position',[315,32,365,48]) 


add block('built-in/To 
Workspace',[sys,'/',['Reactivity',13, Model/Drum',13, 'Reacti 
vity']]) 
set param([sys,'/',['Reactivity',13,'Model/Drum',13, 'Reactiv 
ity I. 

mat- name T ie Dee 

'buffer','40000',... 

'position',[380,177,430,193]) 





add block('built-in/To 
Workspace',[sys,'/',['Reactivity',13,'Model/Net',13,'Reactiv 
ty) 
set param([sys, '/',['Reactivity',13, 'Model/Net',13, 'Reactivi 
villa 

'mat-name','r net',... 

abutter’, 40000 sar. 

'position',[540,25,630,55]) 


add_block('built-in/Note',[sys,'/',['Reactivity',13,'Model/N 
et Core Reactivity']]) 
set_param([sys,'/',['Reactivity',13,'Model/Net Core 
Reactivity']],... 

gposition',[550,175,555w180]) 


add _block('built-in/Outport',[sys,'/',['Reactivity',13,'Mode 

1/Net dk']]) 

set_param([sys,'/',['Reactivity',13,'Model/Net dk']],... 
'position',[605,250,625,270]) 


add block('built-in/Saturation',[sys, /',['Reactivity',13,' M 


odel/Saturation']]) 
set param([sys,'/',['Reactivity',13, 'Model/Saturation']],... 
over Limit','-5!,... 

QUpper Limit' 'Beg*.9',... 

'position',[525,248,555,272]) 
add_line([sys,'/',['Reactivity',13,'Model']],[265,95;310,95; 
310,130;360,130]) 
add_line([sys,'/',['Reactivity',13,'Model']],[295,185;310,18 
5;310,140;360,140]) 
add line([sys,'/',['Reactivity' ,13, 'Model']],[70,95;225,95] ) 
add _ line([sys,'/',['Reactivity' ,13, 'Model']],[280,95;280, 40; 
310,40] ) 
add_line([sys,'/',['Reactivity' ,13, 'Model']],[295,185;375,18 
5] ) 
add line([sys,'/' 
5;455,40;535,40]) 
add line([sys,'/' 
]) 
add_line([sys,'/',['Reactivity',13,'Model']],[560,260;600,26 
0]) 
add_line([sys,'/',['Reactivity' ,13, 'Model']],[455,135;465, 26 
0;520,260]) 


,['Reactivity',13,'Model']],[390,135;455,13 


,['Reactivity',13,'Model']],[65,185;200,185 
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% Finished composite block ['Reactivity',13,'Model']. 


Se 13, Modelli... 
'position',[190,69,220,116]) 

add line(sys,[135, 145; 140, 145;140,105;185 ‚105]) 

add_line(sys,[360,95; 535 ‚95]) 

add_line(sys,[415,95;415,255; 440,255] ) 

add line(sys,[415,155;435,155]) 

add_line(sys,[575,95;585,95;585,50;165,50;165,80;185,80]) 

add_line(sys,[225,95;245,95;245,185;275,185]) 

add_line(sys,[235,95;320,95]) 

add_line(sys,[416,255;40,255;40,145;95,145]) 


drawnow 


% Return any arguments. 
if (nargin | nargout) 
% Must use feval here to access system in memory 
if (nargin > 3) 
if (flag == 0) 
eval(['[ret,x0,str,ts,xts]=',sys,'(t,x,u,flag);']) 
else 
eval(['ret =', sys,'(t,x,u,flag);']) 
end 
else 


[ret,xO,str,ts,xts] = feval(sys); 
end 
else 
drawnow % Flash up the model and execute load callback 
end 
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